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Abstract 

This report presents the performance of a steady-state, dual-spool, separate-exhaust turbofan engine, 
with an inter-stage turbine burner {ITB) serving as a secondary combustor. The ITB, which is located in 
the transition duct between the high- and the low-pressure turbines, is a relatively new concept for 
increasing specific thrust and lowering pollutant emissions in modem jet-engine propulsion. A detailed 
off-design performance analysis of ITB engines is written in Microsoft® Excel macrocode with Visual 
Basic Application to calculate engine performances over the entire operating envelope. Several design- 
point engine cases are pre-selected using a parametric cycle- analysis code developed previously in 
Microsoft” Excel, for off-design analysis. The off-design code calculates engine performances (i.e., thrust 
and thrust-specific-fiiel-consumption) at various flight conditions and throttle settings. 


Nomenclature 

A Flow area 

a Sound speed 

C Constants 

c p Specific heat at constant pressure 

F Uninstalled thrust 

/ Fuel/air ratio 

g c Newton’s constant 

go Acceleration of gravity 

h Altitude 

hpR Low heating value of fuel 

h t Total enthalpy 

K Constants 

M Mach number 

m Mass flow rate 

N Engine speed 
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h m 

J 
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cL 
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prod 
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Rspec 
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std 

TH 

t 

tH 

tL 


CSH 

FTP 


Total pressure 

Gas constant or radius 

Uninstalled specific fuel consumption 

Temperature 

Total temperature 

Absolute velocity 

Molecular weight 

Geo-potential altitude 

Bypass ratio 

Ratio of tree-stream total temperature to SLS temperature of standard day 
Efficiency 

Mechanical Efficiency 
Ratio of specific heats 
Total pressure ratio 
Total temperature ratio 
Total enthalpy ratio 


Subscripts 


Main burner 

Properties between upstream and main burner, engine core or compressor 

High-pressure compressor 

Low-pressure compressor 

Diffuser 

Fan 

High-pressure spool 
Layer of atmosphere 

Properties between exit ITB and downstream, or ITB 

Low-pressure spool 

Maximum 

Mixture of combustion product with air 

Nozzle 

Fan nozzle 

Overall 

Propulsive 

Combustion product 

Reference values 

Military specification 

Ram 

standard properties 
Thermal 

Properties between main burner exit and ITB, or total properties 
High-pressure turbine 
Low-pressure turbine 


Abbreviations 


Constant Specific Heat 
Full Throttle Performance 
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HPC High-pressure Compressor 

HPT High-pressure Turbine 

ITB Interstage Turbine Burner 

LP Lapse Rate 

LPC Low-pressure Compressor 

LPT Low-pressure Turbine 

MFP Mass Flow Parameter 

MSH Modified Specific Heat 

PTP Partial Throttle Performance 

SLS Sea Level Static 

TR Throttle Ratio 


1. Introduction 

This engine performance analysis relies heavily on and builds upon the engine performance analysis 
in references 1, 2, and 3. As eloquently stated in reference 2, the parametric cycle analysis of the previous 
study (refs. 4 and 5) showed “how reference point ( design point or on-design) engine performance was 
determined by design choices.” This made it possible to examine trends in engine specific performance 
{specific thrust and thrust specific fuel consumption) with changes in design variables and to begin to 
narrow the desirable range for each design parameter. Now that the first step has been taken toward 
finding an optimum engine for a particular application, it is time to determine an engine’s steady state 
operating performance. 

The objective of performance cycle analysis is to determine the engine’s performance over its 
operating envelope. The performance of several different engines can then be compared to reveal trends 
in engine performance with design choices. It is thus possible to focus on the most promising design and 
to ultimately find the engine reference point that has the best-balanced performance over the entire 
mission spectrum. 

In general, performance analysis differs significantly from parametric analysis. In parametric 
{reference or on-design point) cycle analysis, all of the design choices (including the flight conditions) are 
free to be selected by the designer, and the specific engine performance characteristics are determined for 
each selected set of choices. In contrast, in performance cycle analysis or off-design analysis, the design 
choices have been set, and the performance of this specific reference point engine is then analyzed at all 
possible operating conditions. The independent variables in performance analysis are different flight 
conditions, throttle settings, and nozzle settings. Once the engine is sized and the mass flow rate is 
specified, performance analysis is employed to determine how a selected engine performs at all operating 
conditions within its flight envelope. 

At the stage of performance analysis, it is considered that an engine of fixed configuration has been 
manufactured for testing in a facility that can simulate the desired flight conditions. Although the absolute 
size of the engine may be arbitrary during performance analysis, and we may lack most of the details of 
the geometry of each individual component, it is still best to think of it as pieces of hardware that have 
been built in order to explore the behavior of this reference point engine. 

In off-design performance analysis, there are two classes of predicting individual component 
performance. First, each individual component performance can be modeled as a function of operating 
conditions. Secondly, actual component characteristics can be obtained from component hardware 
performance data for better estimate of the expected engine performance. These two approaches 
correspond to varying degrees of accuracy. In the absence of actual component hardware in a preliminary 
engine design, simple models of component performance in terms of operating conditions are used to give 
preliminary estimate of the engine performance. 
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In the following section, relationships are developed for determining the engine performance of the 
engine cycle whose parametric ( reference point ) performance was covered in the previous study (refs. 4 
and 5). The analysis is an extension of that presented in references 1, 2, and 3. The task now is to develop 
design tools in the form of equations and their solution procedures that will permit determination of the 
component, and hence engine, performance at any condition. 


2. Aircraft Engine Performance Parameters 

2.1 Performance Parameters 

2.1.1 Thrust . — Thrust is the force used to sustain the flight (thrust = drag), accelerated flight (thrust > 
drag), deceleration (thrust < drag). Using figure 3.1 for the control volume, we can apply momentum 
balance to the control volume. Uninstalled thrust F of a jet engine (single inlet and single exhaust) is 
given by 




Sc 


(la) 


where m {) , m j uel are mass flow rates of inlet air and fuel respectively 
V 0 , V e are velocities at inlet and exit respectively 
P 0 , P e are pressure at inlet and exit respectively 

2.1.2 Specific Fuel Consumption . — It is the rate of fuel use by the propulsion system per unit of 
thrust produced. The uninstalled specific fuel consumptions, S, are given by 

^ ™ fuel 

F 


2.2 Notations 


Some useful quantity notation for compressible flow will be used in this report namely, stagnation 
temperature, stagnation pressure, and Mach number. 

Stagnation temperature or total temperature T t is defined as the temperature obtained when steadily 
flowing fluid is brought to rest adiabatically without involvement of works. Applying first law of 
thermodynamic to calorically perfect gas gives: 


where, 


K = h t = h t 



h t = enthalpy at stagnation condition 

With the assumption of constant specific heat coefficient, the above equation can be written as: 


T t - T + 


V 

2 C, 


T, = T 


1 + 


Y-l 


M' 
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where, 


T t 

T 

M 

J 


stagnation temperature 
static temperature 
Mach number 
ratio of specific heats 


Stagnation pressure or total pressure Pj is defined as pressure reached when a steady flowing fluid is 
brought to rest adiabatically and reversibly. Using the isentropic relation, total pressure is given by 


f 


P t =P 


1 + 

V 2 


Y 1 M 2 


W-i 


Ratio of total temperature (x) and ratio of total pressure (7x) across a component is denoted by 
subscript: d for diffuser, cL for low pressure compressor, cPl for high pressure compressor, b for main 
burner, itb for inter-stage turbine burner, tL for low pressure turbine, tH for high pressure turbine, n for 
nozzle, and / for fan. 

For example: 

total pressure leaving high - pressure compressor 

n cH — 

total pressure entering high - pressure compressor 

total temperature leaving high - pressure compressor 

x cH — — 

total temperature entering high - pressure compressor 


For the free stream, ram, we define x r as a ratio of total temperature/static temperature and n r as a 
ratio of total pressure/static pressure. 


x„ = ^_ = 1 + ^ — -M 


Y-l 


no 


1 + — — - Afn 


y/(y-1) 


A ratio between total enthalpy of the burner exit and ambient enthalpy, denoted by x%, is defined such 
that it will be one of the input parameters. 


x l-b 


(c c„T t ) 

v y ' burner exit 

[c D T t ) ~ 
v y 1 ' ambient 

Ur), 


x \-itb 


"P > ’ITB exit 




P * ' ambient 
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2.3 Component Performance 


In this analysis it is acceptable to assume that the working fluid in the engine can be idealized as a 
perfect gas. Properties of an ideal gas strongly depend on the temperature. This cycle allows fluid 
properties variation across the engine which assumes constant fluid properties from the main burner 
entrance upstream (c pc , y c ), from ITB entrance to main burner exit (c pt , y t ), and from ITB exit downstream 
(Cpt, Yt)- 

2.3.1 Inlet and Diffuser. — Pressure losses occur due to the friction with the inlet wall. The pressure 
total ratio, kj, is always less than 1. 

In supersonic flight, the pressure losses cause shock waves which produce greater pressure losses. 

The inlet total pressure is defined as the product of the ram pressure ratio and the diffuser pressure ratio. 
Therefore the portion of the pressure losses due to the shock waves and wall friction is defined by: 

^ d ~ ^ d max dr 

From the Military Specification 5008B (ref. 6), the following relation is obtained: 

f 1 for Mq < 1 

T[R |i-0.075(m o -l) 1 ' 35 for \<M 0 <5 


2.4 Parametric Cycle Analysis 

The performance of a selected reference point of a turbofan engine is desired at any flight conditions, 
throttle settings, and nozzle settings. It is assumed that a parametric ( design point) cycle analysis has been 
performed for the reference point engine using the method of the previous study (ref. 4) to give the so- 
called reference conditions (subscript R) for the engine (Sr, [F/m 0 \ R , etc.), for each engine component 

(kjr, XfR, etc.), and for the flight conditions (Mqr, Pqr, and Tor). 

To better understand the methodology of studying the performance cycle analysis problem, it is 
instructive to review how the parametric analysis proceeds. For example, in the parametric cycle analysis 
problem as provided in (ref. 4), there are 1 1 independent equations for solving the values of the 1 1 
dependent component parametric variables (e.g., y, x c l, x r /[,f„ x t i, n t H, x t i, n t L,ftb, Mi>, and M\f). These 
1 1 equations are solved for the given values of the independent quantities consisting of flight conditions 
(i.e., Mq and 7o) and design choices (i.e., 7 if, k c r, n c n, and a). With values of these 1 1 component 
variables in hand, the engine reference point performance in terms of F/m 0 , S, r\p, and v\th are readily 

found. A more complete description of the on-design code is given in reference 4 which is published with 
this manual. 

Similarly, to find the off-design engine performance, the operational performance values of the 1 8 
dependent variables listed in table 1 must be determined and then, in turn, related by 1 8 independent 
equations to obtain a solution for each component performance variable. Since performance analysis is an 
indirect problem as opposed to the direct problem of parametric analysis, the solution of the 1 8 
performance equations is not as straightforward as in the parametric case. Regardless of the difficulty in 
solving the equations, once the values of the 1 8 dependent variables in table 1 are known, the engine 
performance in terms of F / m 0 , S, r|p, and v\th follows immediately. In the “Constant or Known” 

column of table 1, 7i/„ x it /,, n n f, and n lt are assumed to remain constant, and xj, x n , and x„yare assumed 
equal to one. 
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Table 2.1. — Engine Performance Variables 


Component 

Independent 

Constant or Known 

Dependent 

Engine 

Diffuser 

Mo, To, Pq 

max o) 

m 0 , a 

Fan 


n/ 

n i x i 

Low-pressure compressor 


TlcI 

TlcL-> cL 

High-pressure compressor 


Tied 

KcH, x cH 

Main burner 

T , 4 

Tib 

fb 

High-pressure turbine 


AtH 

KfH, x tH 

Inter-stage turbine burner 

Tf 4.5 

Kith 

fitb 

Low-pressure turbine 


T ML, A5 = f %itb T 

A.5 R \*itbR J 

KtL, x tL 

Core nozzle 


n„, -^ = f— T 
A/! V X itbR J 

Mg, Po/Pg 

Fan nozzle 


Tl n f 

Mig, Po / Pl9 

Total number 

5 


18 


3. Separate Exhaust Turbofan with ITB Cycle Analysis 

3.1 Assumptions 

The type of engine considered in this study is a separate-exhaust, two-spool turbofan engine with an 
Interstage Turbine Burner {ITB) serving as a secondary combustor as shown in figure 3.1. ITB, is located 
at the transition duct between the high- and the low-pressure turbines, is a relatively new concept for 
improving engine performance. 

The station numbers of each location as indicated in figure 3.1 are in accordance with Aerospace 
Recommended Practice (ARP) 755A (ref. 6). 



2.5 3 4 4.4 4.5 5 7 8 9 


0 2 13 1819 

Figure 3.1 . — Station Numbering for Separate-Exhaust Turbofan Engine 
with Inter-stage Turbine Burning (ITB) (ref. 6) 
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In addition to the assumptions summarized in the parametric analysis (ref. 4), the following 
assumptions are employed in performance analysis: 

• The flow area is constant at station 4, i.e., main burner exit or high-pressure turbine entrance; 

• The flow is choked at the high-pressure turbine entrance nozzles (station 4), at the low-pressure 
turbine entrance nozzles (station 4.5), and at the throat of the exhaust nozzles (stations 8 and 18). 
Since the exhaust nozzles may unchoke at low throttle settings and influence the fan operating 
line, the cases of the unchoked exhaust nozzles (stations 8 and 18) are also included in this 
analysis; 

• The component efficiencies (r\j, r\ cL , r\ cH , r\tL, Tz>, i Vtb, r imfh and r\ mL ) are constant; 

• The total pressure ratios (n ( / max , n/,, nm„ tc„, and k„j) remains the same as the values in (ref. 4). 

• The simplifying gas model of a calorically perfect gas is used in the analysis. It assumes that the 
gases are calorically perfect upstream and downstream of the burner and inter-stage turbine 
burner (ITB) and values of y t , c t , Jm,, and cm, do not vary with throttle setting. For this model, 
variations in fuel-air ratios (f, and f t b) are ignored when compared to unity. 

• There is no turbine cooling. 

• The exit areas (Ag and A \gj of the exhaust nozzles are adjustable so that the corresponding 
pressure ratios ( Pg/Pg and P\g/Pg) can be set to predetermined values. 

• The area at each engine station is constant. However the areas of stations 4.5 and 8 change with 
the ITB setting to follow a specified function of xjtb- 

• The diffuser total pressure ratio, Km is given by 


Kd max Ti? 


where Kd max is the total pressure ratio due to friction and rj R is the total pressure recovery for the 
shocking system given by MIL-E-5008B (ref. 7). 


3.2 Referencing 

The functional relations for engine cycle analysis are based on the application of mass, energy, 
momentum, and entropy considerations to the one-dimensional steady flow of a perfect gas at an engine 
reference or off-design steady state operating point. Thus, if, at any off-design point, 

f{x,Ti) ~ constant (3.1a) 

represents a relationship between the two performance variables x and 7t at a steady state operating point, 
then the constant can be evaluated at the reference point, so that 

f(x,n) = f(x R ,n R )= constant (3.1b) 

sincc /(x,Ti) applies to both on-design and off-design points. This technique for replacing constants with 
reference point values is used frequently in the performance analysis to follow. For example, see eq. (3.9). 
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3.3 Mass Flow Parameter 


The mass flow parameter ( MFP ) is derived by combining mass flow per unit area with the perfect gas 
law, Mach number, sound speed, and equations for total temperature and pressure. The resulting 
expression is 


MFP(M) = 


P t A 



(3.2) 


The total pressure mass flow parameter may be used to find any single flow quantity when the other 
four quantities and the calorically perfect gas constants (y and R) are known at that station. The MFP is 
often used, for example, to determine the flow area required to choke a given flow (i.e., at M= 1). The 
MFP can also be used to develop valuable relationships between the flow properties at two different 
stations, especially when the mass flow is conserved between them. 


3.4 High-Pressure Turbine 


For our engine model, the mass flow entering the high-pressure turbine equals that entering the main 
burner plus the fuel added in the main burner. Thus 

m 4 = m 3 + rii h = /h 3 ( 1 + f b ) (3.3) 


Likewise, the mass flow entering the low-pressure turbine equals that entering the high-pressure turbine 
plus the fuel added in the ITB. Thus 

m 45 = m 4 + m itb = m 3 (l + f b + f itb ) (3 .4) 

Writing these using the mass flow parameter and the flow properties at stations 4 and 4.5 gives 

m 4 =^-A 4 MFP{M 4 )= »’? 3 (l + f b ) (3.5) 

1 tA 


and 


«4. 5 = ^ A 4 5 MFP(M 4 5 ) - 777 3 (l + f b + f itb ) (3.6) 

2/4.5 

Solving each of these two expressions for 777 3 and equating yields 

P t4 A 4 MFP{M 4 ) _ P t4 5 A 4 5 MFP(M 4 5 ) 

2/4 (l + fb ) 2)4.5 + fb + fitb ) 

which can be rewritten in terms of ratios as follows 

P t4A J _ A MFP(M 4 ) (1 + fb + fitb ) P t \A p g, 

2)4 45 4 MFP{M 45) (1 + f b ) P t 4.5 

The right-hand side of the above equation is considered a constant since we make the following 
assumption: the flow is choked at stations 4 and 4.5, the flow areas at station 4 is constant, variations in 
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fuel-air ratios (//, and /,7/,) are ignored when compared to unity, and that the total pressure ratio of the ITB 
does not change. Using the referencing, we have 


r t4.4 




14 


Jr, 


A 4.5 ~ 


74.5 


1 4.4 


-jTt4 


n4 




1 4.5 


? 4.5 


(3.9) 


which can be rewritten in terms of component total pressure and total temperature ratios as 

■\l t tH~t it b A 45 R 


KtH 


x tH K itb , f> ^4.: 


~n t HR 


(3.10) 


This gives one relationship between the high-pressure turbine’s total pressure and total temperature ratios 
to meet mass conservation as a function of the ITB ' s total temperature ratio. A second equation between 
Kin an d Kin comes from the high-pressure turbine efficiency equation or 


*tH 1/Y ' 


(3.11) 


For reference values of n t jj, K t n, A 45 , and %/,, and the corresponding value of v\th, equations (3.10) 
and (3.11) yield the high-pressure turbine’s operating point (r,//, it/#) for a given value of i/,/, and /I 4 5 as 
sketched figure 3.2. The operating point is the intersection of equations (3.10) and (3.11). 

One relationship for A 4 . 5 /A 4 . 5 R is to assume that it is related to the total temperature ratio of the ITB 
(Kith) raised to the power n or 


^4.5 /^4.5 R 



(3.12) 


For the case where A 4 5 /A 43 R - 
n tH = n tHR an d K tH = x tHR . 


li^ith / K itb R ) 


, then the solution to equations (3.10) and (3.1 1) is 



Figure 3.2. — Variation in Fligh-Pressure Turbine Operating Conditions with Increasing ^ x itb / A 4 < 
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3.5 Low-pressure Turbine 


The mass flow rates at station 4.5 equals that at station 8 . Writing the mass conservation at stations 
4.5 and 8 using the mass flow parameter and gives 


-^a 45 mfp{m 45 ) 

VT'4.5 



a 8 mfp(m 8 ) 


(3.13) 


This can be rewritten as 


71 IL _ A 4.5 MFP(M 4s ) 

n n A 8 MFP(M 8 ) 


Noting that n„, and M 4 5 (=1) are assumed constant, then the above equation can be rewritten using 
referencing as 


x lL A 8 R A 4.5 M FP{ M 8 r) 
V x tLR A 8 A 4.5R MFP(M 8 ) 


(3.15) 


This gives one relationship between the low-pressure turbine’s total pressure and total temperature ratios 
to meet mass conservation as a function of the core exhaust nozzle throat area ( A 8 /A 8 r ). A second 
equation between n t L and x t i Comes from the low-pressure turbine efficiency equation or 

x tL =\-^ L {\-n^ L itb ~^ hub } (3.16) 

The operating condition (x t i, n t £) of the low-pressure turbine can be determined using equations 
(3.15) and (3.16) along with the exhaust Mach number Mr and specified functional relationships for the 
low pressure turbine nozzle throat area (^4 5/^4 5 ^) and the core exhaust nozzle throat area (A 8 /A 8 r). One 
relationship for Ag/AgR is to assume that it is related to the total temperature ratio of the ITB (Jnh) raised to 
the power m or 


A 8 /- A 8 R ~ ( T itb / 1 itb R Y 


(3.17) 


With this functional relationship, the engine’s low-pressure turbine performance will vary the same as 
the turbofan without the ITB, when the ITB is turned off. For reference values of n t i, x t L, and x,,/, and the 
corresponding value of r\ t i, equations (3.15) and (3.16) yield the low-pressure turbine’s operating point 


( x t L , n t L ) for a given value of 


A 8R A 4.5 MFP(m 8R ) 


a 8 a 4 . 5 r MFP{m 8 ) 

intersection of equations (3.15) and (3.16). For the case where A 8 /A 8 r 
solution to equations (3.15) and (3.16) is n t i = n t LR and x t i = x t if>. 


as sketched in figure 3.3. The operating point is the 
A 4 ^/A 4 s r and M 8 /M 8 r, then the 


a 8 a 4 5R mfp{m 8 ) 

A 8R A 4.5 MFP(M 8 r ) 
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Eq. (3.16) 



Eq. (3.15) 


Increasing 


A Am< MFP(M h ) 
Ar 4 *. s mfp{Mz R ) 


Figure 3.3. — Variation in Low-Pressure Turbine Operating Conditions with Increasing 


3.6 Engine Bypass Ratio 

An expression for the engine bypass ratio (a) follows directly from its definition written in terms of 
properties at stations 4 and 18. 


where 


4 T tl% 


-MFP l8 = 


P 0 n r n d n j-n n fA 18 
V^13 


m 4 = m c (l + f b )= *aA^MFP 4 = P O n r^cL^cH^bA MF p 4 
yJ P t4 V^4 


Thus the engine bypass ratio can be written as 


(1 + A) f^7 MFP{M x 8 ) A 18 

^cL^cH^b4 1t f 1t nf)i r A3 MFP{M 4) A 4 


Since ( I +//,), 7i/„ n n f,A\%, A 4 , and M 4 (=1) are assumed constant and equal to their reference values, the 
bypass ratio can be rewritten using referencing as 


a - a r 


n cLR n cLR/ n fR T t 4 /T t4R MFP(m 18 ) 
n cL n cH l n f Kv/( T r* T /R) MFP ( M 18 «) 
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3.7 Fan and Low-pressure Compressor 


The performance equation for the total temperature ratio of the fan follows directly from the power 
balance of the low-pressure spool. Equating the power received by the fan and low-pressure (LP) 
compressor to that delivered by the low-pressure turbine gives 

™ f c ~ Tf2 ) T W-cCpc (7;2.5 — ^t2 ) — r \mL^4.5 c pitb{^t4.5 ~ (3.23) 

or 

^f c pc^t2 { a (j/ ~^)~ >r ^LcL — l)j — T fb f itb^pitb^t 4.5^ ~^tL) (3.24) 

Then dividing by m c c pc T t2 , gives 

a{xf - ])+ (x cL — l) = ri mL^ + fb + fitb ) X Ub 0 ~ x il) (3.25) 

T,. 


The temperature rise of the low-pressure compressor can be approximated as proportional to that of the 
fan (refs. 2 and 3) and thus 


let 

V 


-1 


*cL 

Tf-l 


(3.26) 


Solving equation (3.26) for (x c p - 1), substituting into equation (3.25), and then solving for the fan 
total temperature ratio (xj), gives 


= 1 + (Vr “ ! )l 


1 mL ' 


T X- - itb_ J ~ x tL Xj_ + fb + fitb )_ 

“ r 1 *cLR ~ 1+ a(t fr - l) J 

Using the fan efficiency, the fan total pressure ratio is given by 

nf - {l + x\f (t/-i)}^/^ c ^ 

The low-pressure compressor total temperature ratio follows from equation (3.26) and thus 


XcL 


= 1 


'' r til 


U 1 ( 

Using the LP compressor efficiency, the LP compressor total pressure ratio is expressed as 

= {l + hcZ, (xcZ, _ l)} Yc// ^ c ^ 


(3.27) 


(3.28) 


(3.29) 


KcL : 


(3.30) 


3.8 High-Pressure Compressor 

The performance equation for the total temperature ratio of the high-pressure (HP) compressor 
follows directly from the power balance of the high-pressure spool. Solving equation (37) of the 
parametric analysis (ref. 4) for the total temperature ratio of the high-pressure compressor (x c h) gives 


NASA/TM— 2005-213659 


13 



x cH ~ l + r \mH^ + f b) 


T A,-fe(l~Thff) 


(3.31) 


X r X cL 


Using the HP compressor efficiency, the HP compressor total pressure ratio is given by 

^//-{l + hc//( x c //-l)} Yc/(Yc ' l) (3-32) 


3.9 Exhaust Nozzles 

The total to static pressure ratio of the core exhaust nozzle is given by equation (3.31) of reference 4 
as 


P '9 _ P ') 

— ~— n r n d n cL n cH K b K tH K itb n tL n n 
>9 ‘9 

The exit Mach number follows directly using 


(3.33) 


Mg = 


Y itb~ l 


Y ub 


-l 


£ 79 

\ P 9 J 


Yub 


-1 


If Mg > 1 , then M% = 1 , else M% = Mg. 


(3.34) 


The equation for the area ratio of the core exhaust nozzle follows from mass conservation between 
stations 8 and 9 and use of the mass flow parameter, giving 

A g 1 MFP(M 8 ) 

A 8 k„ MFP(Mg) ( ’ 

The total to static pressure ratio of the fan exhaust nozzle is given by equation (3.36) of reference 4 as 

P t\9 _ p 0 „ „ „ „ j i i 

~A~-~A~ n r n d n f n nf (3.36) 

M9 M9 

The exit Mach number follows directly using 


M \ 9 - 


I Y c -1 


-l 


f P 1 

' 119 

V ^19 J 


-1 


\fM\g > 1, thenMi8 = 1, elseMjs = M\g. 


(3.37) 


The equation for the area ratio of the core exhaust nozzle follows from mass conservation between 
stations 1 8 and 1 9 and use of the mass flow parameter, giving 
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(3.38) 


A ig _ 1 MFPjM ig) 

A 18 n nf MFP{M l9 ) 

3.10 Engine Mass Flow Rate 

An expression for the overall engine mass flow rate ( m 0 ) follows by the use of the mass flow 
parameter at station 4, giving 


m 0 - m 0 


1 + q P 0 K l .n d n cL n cH T t 4R 
1 + a^ {Pq K r T^ d ^ cL n cH )r V ^4 


(3.39) 


3.11 Combustor Fuel-to-Air Ratios 

The fuel-to-air ratios for the main burner (//,) and the ITB (f, b ) are determined by the user’s choice. 
There are two models of the gas properties available in this analysis, i.e., Constant Specific Pleat ( CSH) 
model and the Modified Specific Heat ( MSH) model. Both models include the products of combustion in 
the heat capacity. For the details, see appendix B. 


4. Summary of Equations 

This section presents the complete performance cycle analysis equations for non-afterburning, 
separate exhaust flow, two-spool turbofan engine with inter-stage turbine burning. 

Inputs 


Performance choices 

Flight parameters'. Mq, h [T 0 = Tamhf) (K, °R), P 0 =P am b(h) (Pa, psia)] 

Throttle settings : T t 4 , T t 4 5 (K, °R) 

Engine control'. TR 

Design Constants 

Perfect gas constants'. y c , y t , y itb , c pc , c pt , c pit b [kJ/(kg-K), Btu/(lbm-°R)] 

Tt S. Kd maxi ti b , Tt jf b , ti n , tl n f 

H t\f> hcT.’ 9 9///, 9 tL> 9/>’ 9(76’ 9«i//’ 9 mL 
Fuel: hpR (kJ/kg, Btu/lbm) 

Reference conditions 

Flight parameters: M 0R , h R (km, kft)[T 0R = T amh (h R ) (K,°R), P 0R =P amb (h R ) (Pa, psia)], x R , n R 

Component behavior: nj R , n c lr, ti c HR, n t HR, n t LR, fR, tcLR, t cHR , t-tHR, t-tLR 

Others: T t4R , T t45R (K, °R), f bR , f tbR , M 8R , M 18R , a R , F R (N, lbf), m 0R (kg/sec, lbm/sec), S R 

/ mg / s Ibm / hr ' 

ibf ; 
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Engine Control Limits'. 


^(4 max’ ^(4.5 max' tZ f max ’ max >^3 max (^> ^t3 max psia), /oNpp, /oNfjp 

LP Turbine Nozzle Throat Area Functional Relationship'. 

A4.5 / A A 5R = ( x itb i x itbR )" 

Core Exhaust Throat Area Functional Relationship-. 

^8 ! A 8 R = ( tith / litbR 


Outputs 


Overall Performance 

F (N, lbf), m {) (kg/sec, lbm/sec), S 

T9/T0, 


f mg / sec Ibm / hr A 


N 


lbf 


, fo , b p, b THi ba VfaQ, Vifao, a, P&/P9, 


Pt\9lP\% T\ 9 /Tq 

Component Behavior 

Kf, n cL , n cH , n tH , n tL , Xf, x cL , x cH , x tH , x tL , x x-b, ^x-ith ffr fth M 8 , M \ 8 , M 9 , M X9 , A 9 /A 8 , A 19 /A 18 


Reference Conditions 
If (SI unit ) then 


Else 


Equations 


T ref = 288.6 K 


End If 


T ref = 516.67 °R 


X rR ~ 1 + ~~ Mofi 


-.Yc/itc- 1 ) 


^ rR x rR 


*1 RspecR = 1 for Mq R < 1 

^RspecR - 1 _ 0.075(m or -l) 1 " ~~ for 1<M 0 ^<5 


(4.1R) 

(4.2R) 

(4.3R) 
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ft dR ft d max 'H RspecR 


(4.4R) 


ft itbR 


1 I4.5R 

TfARfttHR 


Tift+1 


MFP(M 8R ) = M 8R J yitbSc 


1 'itb \y , ^2 ^2(l y,' ( £,) 


R; 


itb 


1 + - 


8 R 


mfp{m 18R )=m 18R ^ 

V 


y c +i 


' Yc — 1..2 >(1-Tc 


l + 4£ M- 


18 fi 


A __ 

O 0« - T ri? — 

1 ref 


If 0 Oi j > 77? then 


Else 


Endif 

Preliminary computations 1 


Y?4max 04/? 
^/4.5max — ^tA.SR 


04 max 


04.5 max 


-Y?4i? 

V 


- Tf4.5R 


9q R ' 
TR , 

f 9(XR 


A 


y 


x,. = 1 + ^M n 2 
2 0 

n- _ T Yc/(y c -l) 

l\ty — y 


07 ? = 1 M Mq < 1 

r|yj = 1 - 0.075(M 0 -1) 135 for 1<M 0 <5 


(4.5R) 

(4.6R) 

(4.7R) 

(4.8R) 


(4.9R) 
(4.1 OR) 

(4.1) 

(4.2) 

(4.3) 


ftd — ftdvtmf\R 


(4.4) 


R 


c 


Yc- 1 


C 


pc 


(4.5) 
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(4.6) 


Ritb= lM ^-C pt (4.7) 

Y ub 

«o =ylr c Agd D (4.8) 

^0=^ c/(Ye " l} (4-9) 

T t0 =T 0 x, (4.10) 


Control limits 


mfp(m 0 )=m 0 


YcSc 


R, 


1 + 



M 0 2 l 


Tc+l 

2(1-Yc) 




9o 


= V 



If 0 q > TR then 


T = T 

± t 4 ± t 4 max 


Else 


T,a 


- T, 


1 4 max 


TR 


\ 

) 


(4.11) 


(4.12) 


(4.13a) 


(4.13b) 


Endif 


If (A. 1TB on) then 
If 0 q > TR then 


Else 


'/4.5 


= T, 


1 4. 5 max 


2^4.5 ^/4.5max 


0o. 

\TRj 


Endif 

Endif 


(4.14a) 


(4.14b) 
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Set initial values : 


7« 0 - m 0R , f b - f br , f itb - f itbR , 

n tH ~ n tHR > XtH ~ XtHR , KtL = TttLR , XfL = XtLR 

n f ~ n fr> n cL = n cLR’ n cH = n cHR 

x f ~ x fR’ x cL = x cLR > x cH ~ x cHR 


Inter-stage Turbine Burner 

1 

If ( A.ITBon ) then 


x X-b 


CpAt4 

CpAo 


x itb ~ 


x X—irb 


7)4.5 

7)4 x tH 
C pitb T t4.5 

C pc To 


Else 


Xitb = 1 


(4.15) 


(4.16a) 

(4.16b) 


(4.16c) 


_ C p ,T t 4 _ 

XX-itb ~ rp x tH 

'~'pc*0 


(4.16d) 


7 ) 4.5 - Tt4 x tH 


(4.16e) 


Endif 


^4.5 / A.5R ~ (' x itb / x itbR)' 


(4.17) 


A ! A R — ( x itb ( x itbR )” 


High-pressure turbine 


x tH =1 — T] tH ^ /y ' 




sl x tH x itb A.5R 
{^ x tH x itb )* ^ 4.5 




Is 7i 1 i[\j - Tim- < 0.00001? If so, then continue. If not, set n t n - 7 i,hn an d return to 1 


(4.18) 


(4.19) 

(4.20) 
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High-pressure compressor 
2 


x cH = 1 + r \mH ( 1 + fb ) % — 

x r x cL 

n C H = {l + T1 cH (t cH -l^Ar 1 ) 


Main burner 

If CSH model then 


fb — 


^ X-b tj- x cL x cH 
T\bhpR / K ~ x X-b 


Else 


Tf3 - Tq x ,- x cL x cH 


ENTHALPY {T t3 , 0, h t3 ) 


Set initial value of fuel/air ratio at station 4, / 4 / = X ''~ h — T ' x cL x cH 

r\b ll PR ! h 0 ~ x X-b 


3 ENTHALPY {T t4 ,f 4i ,h t 4 ) 


^ _ h t 4 - h t 3 
P/jAff ~h t A 

If I fb ~ A/ I > 0 . 0001 , then - fo and go to 3; else continue. 
End if 

Inter-stage Turbine Burner 
If CW model then 


■/»=(! + /<,) 

^itbhpR / «0 - x /,-(Y/3 


Else 


2,4.4 - T tA x tH 


ENTHALPY ( 7 , 4 . 5 , As;, 6 , 4 . 4 ) 


(4.21) 

(4.22) 

(4.23a) 

(4.23b) 

(4.23c) 

(4.23d) 

(4.24a) 

(4.24b) 
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(4.24c) 


Set initial value of fuel/air ratio at station 4.4, f 4 5/ (i + f b ) Xx ~ ith ^ 

^\itb hpit/ho — ^X-itb 

4 ENTHALPY {T t4 . 5 , f 4 . 5i , h t4 . 5 ) 


/,,„=(!+ h )AiiA±i_ 

^itbhpR - h t 4.5 

If | \f b + f itb - f 4 5 -| > 0.0001, then f 45i - f b + f itb and go to 4; else continue. 
End if 


lf fitb = 0 then T nb = Jt’Cpitb = c pt > and R itb = R t 


Fan and low-pressure compressor pressure ratios 

7 If = {l + r|y(xy — l)} 'V'c /(yc— !) 

KcL ={ 1 + 11 cl(*cL -l)}vd(rc-i) 

Exhaust nozzles 


P t\9 
P\ 9 


p o 

K r K c iK j~K n f 
'19 


Pf9 _ P() 

— TC r Tl d It cL^ cH tt tH ^ i]Vi n 

P) P) 


M Q = 


I litb 


-1 


^ p 'N 

*t9 

K p 9 ) 


litb 


Jub 


-1 


If M 9 > 1, then Mg = 1, else Mg = M 9 


mfp(m 8 )=m% 


Y itbSc 
R itb 


( 

1 + 

V 


Y itb 1 
2 


Ml 


Y/rfc +1 


mfp(m 9 )=m 9 


Y itb Sc 


R 


itb 


Y itb 1 
2 


m 9 2 


YiYi +1 
^2(1-^) 

y 


h 9 _ i mfp(m 8 ) 
a^~^mfp(m 9 ) 


(4.24d) 

(4.34) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4-30) 

(4.31) 

(4-32) 
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M\9 = 


(4.33) 


y c ~ l 


' Pm 

V P\9 


Yc " 1 




If Mi 9 > 1, then Mis = 1, else M\ g = M\ 9 


MFP(Mi 8 ) = M 18 


Ycgc 




1 + 


Yc -1 

2 


"ft 

/ 


Yc +1 

2(1-Tc) 


mfp(m 19 )= m 19 


Ye&c 


fl + 


/?, 


Yc -1 

2 


"l 2 9 


Yc + 1 
>(1-Yc) 


J 19 _ 1 MFP{M n ) 
A 18 n nf MFP(M l9 ) 


Engine bypass ratio 


a-a R 


KcLRKcHR / 
’^■cL't^cH ! tt f 


I T t4 /T t4R MFP(Mi 8 ) 
\x r x f l{x rR Xf R ) MFP(M 18r ) 


Fan and low-pressure compressor temperature ratios 


VjN ~ 



(j_ ~ )(l + //) + fitb ) ] 

mL T r — 1 + Cx(xy^ -l) 


i , ,)(^ 

x cL = 1 + - 1)7 rt— 

Low-pressure turbine 

*tL=l~ntL {\-n^ tb ~^ lyitb ) 


I tfL ^8 R ^ 4.5 MFP(M$r) 

htZR A A A , 5R mfp(m 8 ) 


Is iflj - x f < 0.0001? If so, then continue. If not, set x f = x ^ an return to 2 


(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 


NASA/TM— 2005-213659 


22 



Remainder of Calculations 


Engine mass flow (mo) 


m 0 = riiQ R 


1 + CX PqK i -K ( jK cR 71 c }[ T 4 R 
1 + a R (P \) n r n d n cL n cII )/? V T t4 


(4.42) 


Overall engine performance 


/o - (fb + fitb)/^ + a ) 

(4.43) 

T 9 _ T t 45 x tL /T 0 

To (P t9 / p^fiub-iVia, 

(4.44) 

T l9 _ T r Xf 

To {P tl9 / P l9 )^~ 1 ^ 

(4.45) 

a 9 — ■>/ Y itbR-itb§ cT 9 

(4.46) 

a l9 ~ si lc R cgcT \ 9 

(4.47) 

V 9 — M 9 a 9 

(4.48) 

V \ 9 - M\ 9 a \ 9 

(4.49) 


F_ 

m 0 


a o 

g c (l + a) 


[l + /o (l + a)]— + a — ^-(l + a)M 0 
a 0 a 0 

• + [i + /o (i + a )]^ML > 

L R c V 9 /a () Yc 

, a t 19 /t 0 {i-p 0 /p 19 ) 

V\g/a 0 y c 


(4.50) 


F = 777 Q 


^ F A 


777 0 


(4.51) 


5 = 


/o 

F/w 0 


(4.52) 


2 g c A4o(l + a) 


ilP 


F 

777 o 


ao-|[l + /o(l + a)1 — 

\a 0 y 


A2 


+ a 


V a 0 


(4.53) 


- (l + a )Mx 
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dm - 


[l+/ 0 (l+a)]p 
\ao ; 


l 

i Vl9 ) 

+ a 


\ a o ) 


-(1 + a )M 0 2 


2 S c 


A o 


(1 + a) fohpR 

■no = nranp 

~ P t0 MFP(M 0 ) 


(4.54) 

(4.55) 

(4.56) 


%RPM IP spool 


= 100 I Wv-O 

V[jot4v->)], 


(4.57) 


%RPM HPlpcol = 100 r 0 VT cZ (tc 1) (4.58) 

V Fo'c,-x c i (x c// -1)J^ 

If any of the control limits ( 7t c , 7^ , P (3 , etc.) are exceeded, first reduce 7’, 4 5 until 1,7/, = 1 or reduce 
T t 4 and then go to 1. A Newtonian iteration scheme is used to rapidly converge on the values of 7’, 4 5 and 
T t 4 that meets the most constraining control limit. 


5. User’s Manual 

The Microsoft” Excel program is written in combination between spreadsheet neuron cells, visual 
basic, and macro code. These three combinations provide user-friendly software that compilation and 
preprocessing are no longer necessary. The input is to well-labeled Microsoft ” Excel cells. Most of the 
values are prespecified and an example case is displayed. It is good practice to save this case by another 
name before modifying the spread sheet in order to retrieve the starting case. User obtains result plots 
right away just by clicking some buttons. 

The program mainly comprises of three sections, namely, Input, Full Throttle Performance and 
Partial Throttle Performance. The Microsoft” Excel spreadsheets associated with each section are listed 
in the table 2 as shown below: 


Table 5.1. — List of Microsoft 8 Excel Spreadsheets 


No 

Section 

Function 

Microsoft® Excel Spreadsheet 

1 

Input 

Introduction 

Coverpage, Instruction 

Data Input 

Input 





Control 

Control FTP, SinglePt FTP 

2 

Full Throttle Performance 

Data Output 

D FTP, D FvsMO.D SvsMO, D_picHvsM0, 
D Tt4vsM0, D mdotO, D_pif, D alp, 

D Tt4p5 

3 

Partial Throttle 

Control 

Control PTP 

Performatice 

Data Output 

D PTP 
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5.1 Input 


5.1.1 CoverPage. — The ‘CoverPage’ sheet contains the information about the authors of this 
program. Any questions regarding the program can be addressed to the authors through the email or 
phone. 

5.1.2 Instructions. — First-time users are strongly recommended to read this sheet before running the 
program. In this sheet, you will find details of how to run the program. This sheet also explains several 
assumptions made in the equations so that the users are aware of some cases in the equations that have 
been idealized to simplify the problems. 

5.1.3 Input. — This sheet is where most of the inputs are specified. All cells in ‘dark green’ (^ 1 ) 
color need to be filled in with input values, except those cells in ‘hotpink’ (SI) color, where the built-in 
values are set. The program will check all the inputs value in this sheet to make sure that the inputs are 
specified correctly. Window dialog boxes will pop up and notice the user if there are input values that are 
not specified or zero input values are found. 

Once the design choice of a specific reference-point engine has been made from a parametric cycle 
analysis (previously developed in ref. 4), the input values for the so-called ‘ Reference Conditions’ (i.e., 
the flight condition, throttle settings, engine component, etc.) and the ‘ Design Constants’ are readily 
obtained. By clicking the ‘ Import data from 1TB on-design code ’ button, the code will automatically load 
the reference engine data file at a specified path and filename. All reference engine data will be displayed 
at all ‘burlywood’ cells (IH). Please refer to Appendix F for further explanation and instruction on this. 
Alternatively, users can choose to manually input each input value, which is a time-consuming process. 

There are four combo boxes in the ‘Input’ sheet ( Combo Box is a list box that displays a list of values 
and lets the users select one of the value in the list) namely Unit System, ITB Switch, Fuel Burn Model, 
and Perform Computation in as shown in figure 5.1. You need first to specify the value in combo box 
Unit System, combo box ITB Switch, and combo box Fuel Burn Model before moving on to combo box 
Perform Computation in. 

Combo box Unit System lets you specify the input and output unit system. Currently, the program can 
handle two units systems, namely, English and SI units. Combo box ITB Switch lets user to turn ON or 
OFF the Inter-stage Turbine Burner (ITB). This feature provides a flexibility to choose two types of 
engine and they are the engine with ITB-ON and with ITB-OFF. With this feature, you will be able to see 
how much engine performance you can get with ITB-ON and with ITB-OFF. Combo Box Fuel Burn 
Model lets user to choose two models for the computation of the fuel-air-ratio inside the burners, i.e., 
Constant Specific Heat (CSH) model and the Modified Specific Heat ( MSH) model. For the details, see 
appendix B. 

Once all the inputs values are specified and all combo box values on Unit System, ITB Switch, and 
Fuel Burn Model are specified, user can specify a value in combo box Combo box Perform 
Computation in. This combo box provides three option lists to choose from, and they are: 

1. Full Throttle Performance (FTP): Multiple-point 

2. Full Throttle Performance (FTP): Single-point 

3. Partial Throttle Performance (PTP) 

Once you select either one of them, you will be directed to some new sheets. The macro code 
associated with this combo box will open a relevant control sheet, and some hidden data sheets depending 
on the selection you make. 

The following sections describe the instruction of how to run the program for each different option. 
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Table 5.2. — List of Combo Box in ‘ Input ’ sheet 


No 

Combo Box 

List value 

1 

Unit System 

SI 

English 

2 

ITB Switch 

ON 

OFF 

3 

Fuel Burn Model 

CSH (Constant Specific Heat) 

MSH (Modified Specific Heat) 

4 

Perform Computation in 

Full Throttle Performance (FTP): Multiple-point 

Full Throttle Performance (FTP): Single-point 

Partial Throttle Performance (PTP) 


INPUT SHEET 


Reference conditions 


Import data from ITB on-desiqn code | 

Fliaht parameters 


Mqr - 

0.01 


Altitude, hR = 

0.00 

km 

Tor = 

288.15 

K 

Pqr = 

101325.00 

Pa 

Throttle settinas 


Tt4R - 

1450 

K 

Tt4.5R = 

1350 

K 

Component behaviors 


KtHR = 

0.3145165 


*tLR = 

0.602753 


71 cHR = 

8.23 


71 cLR = 

2.43 


*fR = 

2.43 


TtHR = 

0.7899483 


TtLR = 

0.894904 


x cHR = 

1.9434872 


IcLR = 

1.3377882 


TfR = 

1.3144672 


Fuel-air ratios 

CSH 


f bR = 

0,0265885 


f itbR = 

Others 

0.0043636 


MgR = 

1 


MigR = 

1 


O.R = 

0.73 


Fr = 

80962.511 

N 

Sr = 

26.076028 

(mg/s)/N 

mdotoR = 

118 

kg/s 


Performance choices: 

Throttle settings 


T, 4 = 

1450 

Tt4.5 = 

1350 

Exhaust Nozzle 

6 

CD 

Tl 

II 

1.000 

fan: P 0 /Pi9 = 

1.000 

Enaine Control 

JR = 

1.00 

esign constants: 

Total pressure 

ratio 

71 d max = 

0.99 

71 b = 

0.95 

71 itb = 

1.00 

Tin = 

0.99 

71 nf = 

0.98 

Component efficiencv 

ntH = 

0.9112835 

ntL = 

0.9247443 

n ch = 

0.878489 

n cl = 

0.8571859 

nt = 

0.9207551 


Burner efficiency 


rib = 

0.99 

n itb = 

0.99 

Mechanical efficiencv 

nmH = 

0.92 

r|mL = 

0.93 

Low heatina value of fuel 

hpR-b - 

43,124.00 


Cover Page 


Engine control limits: 


Tt4max = 

1450 

K 

Tt4.5max — 

1350 

K 

TIcHmax — 

8.23 


Ttfmax — 

2.43 


Tt3max = 

0 

K 

Pt3max — 

0 

Pa 

%N l = 

100 


%N h = 

100 


Conversion factors: 

9o= 



Conversion = 

1000 


Throat Area Relationship 

n = 

0.5 


m = 

0.5 



Unit system 

fiT 


3 


ITB switch 

I ON 3 

Fuel burn model 

I CSH (Constant Specific Heat) - I 


Remarks: 

I => Input values here. 

=> Do not modify. Data loaded from on-design code 
=> Do not modify. 


Renit 

I 


Reference -point Engine Path: 
h:\reference.txt 


Perform Computation in: 

I Full Throttle Performance (FTP): Multiple-point 


Note : Always check values for Tor and Pqr 

3] FTP Multiple-point 



7 8 9 


02 13 


1819 


Figure 5.1. — Microsoft Excel Input Sheet Screenshot 
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5.2 Full Throttle Performance over Mq and Altitude (FTP) 

5.2.1 Multiple-point Calculations (sheet ‘ControlFTP’). — When this option is chosen, the program 
opens several sheets, namely, one control sheet (i.e., either ‘Control FTP’), and nine data sheets (i.e., 
‘DFTP’, ‘DFvsMO’, ‘D SvsMO’, ‘D_picHvsMO\ ‘D_Tt4vsM0’, ‘DmdotO’, ‘D_pif , ‘Dalp’, and 
‘D_Tt4p5’). 

In ‘Control FTP’ sheet, user expects to see a sheet as shown in figure 5.2. Before the computation, 
user needs to complete two input parameters, i.e., Flight parameters and Altitude, as indicated in ‘dark 
green’ cells. 

NOTE: In this code, there is only one design variable, namely, Flight Mach number. Other iteration 
variables, such as ambient conditions ( Tq and Pq), altitude (/z), total exit temperature of main burner or 
ITB, etc., may be added to the code with minor code modification in the future. 

To specify a range for flight Mach number (Mo), user needs to input the values for the lower limit, 
upper limit, and an increment of Mq. If input lower limit of Mq is zero, the code will internally change it 
from zero to 0.01. 

Click the Altitude (h) combo box to select number of plots for different altitudes. Selected number of 
cells will be created adjacent to th q Altitude (h) combo box. Once all the ‘dark green’ cells are filled, it is 
ready to run the program and generate the plots, simply by clicking Calculate button. Data computed will 
be stored in nine data sheets. 

Clicking the Clear data & plot button will clear up all the computed solutions in data sheets as well 
as the eight plots in ‘Control FTP’ sheet. 


3ff-deslqn Multi-point : Full Throttle Performance over M 0 & Altitude (h) 


© 


Flight parameters 


Mo = | 

0 to 

A Mo = 

0.1 

Altitude 


h=| 


1 3 

Range: 0 

< h < 84km 


h, = |l 0 | km 
h 2 =|l5 |km 


® 


Calculate Clear data & plot 


Computing time: 4 sec 

Instructions III PUT dieet 



Thrust Specific Fuel Consumption Vs MO 



Figure 5.2. — Screenshot of Full-Throttle Performance: Multi-point calculations 
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5.2.2 Single-point Calculations (sheet ‘SinglePtFTP ’). — When this option is chosen, the program 
opens one control sheet (i.e., ‘SinglePt FTP’). This option allows the user to predict the off-design 
performance of the selected engine at a fixed Mq and a specified altitude. As indicated in figure 5.3, user 
just needs to specify a single value for each and altitude. Then, simply clicking the Calculate button 
will run the code. All the input values from ‘Input’ sheet and the solutions will be displayed in the same 
sheet. 


Off-Desiqn Single Point : Full Throttle Performance over MO & Altitude (h) 


M 0 = 
Altitude, h = 
- Ambient: 

Input Constants: 


0 


0 

km 

II 

288.15 

II 

o 

CL 

1 01 325.00 


K 

Pa 


IT B ON 


Calculate 


Instructions 


Clear Data 


INPUT sheet 


Mote : Always check values for T 0 and P 0 


Control Limits: 


Parameter: 


Mach number @ 0, M 0 
Altitude, h 

Temperature @0,T o 
Pressure @ 0, P 0 
Total Temp @4,T U 
Total Temp @4.5,Tti 5 

*I,T, 

*cl 

r 

ItcLi Tel 
A:H. tcH 
%l, TlH 
tttL, XtL 

Control Limit: 

LP Spool RPM (% of design pt) 
HP Spool RPM (% of design pt) 
Mach number @ 8, M s 
Mach number @ 1 8, M, 8 
MB - Fuel/Air Ratio (Q 
ITB - Fuel/Air Ratio (W 
Alpha, a 
Thrust, F 

Spec Fuel Consumption, S 
Mass Flow Rate @ 0, mdoto 
Propulsive Efficiency (%) 
Thermal Efficiency (%) 

Overall Efficiency (%) 


100 

100 

1 

1 

0.026588 

0.004364 

0.73 

80962.51 

26.07603 

118 

0.898424 

34.03478 

0.305777 


^ d mat — 

0.99 

ti) : 

0.95 

it ip- 

1 

Jt,= 

0.99 

Hit= 

0.98 

hb = 

0.99 

n ii) = 

0.99 

hr= 

0.92076 

htH = 

0.911284 

UlL = 

0.92474 

UcH = 

0.87849 

hcL = 

0.85719 

*1 mH — 

0.92 

hinL = 

0.93 

PJP 0 = 

1 

P,JP 0 = 

1 

hpf>-i)= 

43124 KJ/kg 






Tttmax = 

1450 

Tuan®- 

1350 

Tt*na» = 

0 

P t3ma« = 

0 

^cHmat = 

8.230453 

%N h = 

100 

%n l = 

100 





Reference 



Test 


Remarks: 



0.01 



0.01 





0.000 



0 





288.15 



288.15 





101325 



101325 





1450 



1450 





1350 



1348 




1.00007 

1.00002 


1.00007 

1.00002 





0.99 



0.99 




2.43 

1.314467 


2.43003 

1.31447 




2.43 

1.337788 


2.43003 

1.3372 




8.23 

1.943487 


8.23042 

1.94348 




0.314516 

0.789948 


0.31452 

0.78995 




0.602753 

0.894904 


0.60275 

0.89489 




=> Input valu 
=> Do notm 
=> Do notm 


Tt4max/Tt4.5max 

99.9125 

99.9778 

1 

1 

0.02659 

0.0043 

0.73 

80913.2 [N] 

26.0358 [ (mg/s) I N ] 
118.001 [kg/s] 
0.89909 % 

34.0619 % 

0.30625 % 


Figure 5.3. — Screenshot of Full-Throttle Performance: Single-point calculation sheet 
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5.3 Partial Throttle Performance (PTP) 

When this option is chosen, the program will open two sheets, namely, a control sheet (i.e., 
‘ControlPTP’), and a data sheet (i.e., T)_PTP’). 

5.3.1 Control PTP . — In ‘Control PTP’ sheet as shown in figure 5.4, user needs to specify four input 
parameters, i.e.,, Altitude, Min % thrust, Flight parameters, and increment of temperature reduction, as 
indicated in ‘dark green’ cells. 

Click the Flight parameters (Mo) combo box to select number of plots for different Mq. Specified 
number of cells will be created and displayed adjacent to the Flight parameters (Mo) combo box. Once 
all the ‘dark green’ cells are filled, it is ready to run the program and generate the plot. 

Simply pressing the Calculate button, the code will calculate the engine performance from 1 00% 
down to minimum thrust entered in Min % thrust cell. The minimum value of min % thrust is internally 
preset to be 1 0%. If any value lower than 1 0% is found, user will be noticed and asked to re-enter the 
value. 

Other than setting a Min % thrust, user may select to enter a lowest allowable exit temperature of the 
main burner, T t 4 . The default value of 7 ) 4 is 1620 °R (or 900 °K). 


Partial Throttle Performance (PTP) 

© 


Altitude 


h =| 10 | km 


Min % Thrust 


© 

© Flight parameters 




M = 2 




Mi : 

Mo = 


© 


© 


AT=| 20 

K 

Calculate | 

Clear data and plot | 


0.8 

1.2 



Instructions 


INPUT sheet 


Optional : For a converged solution, the lowest allowable T t4 is preset to [ 900 |k 


Partial throttle performance 



- Mach 1 

- Mach 2 


Figure 5.4. — Screenshot of Partial-Throttle Performance 
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5.4 Discussions 


5.4.1 Storing Data . — Whenever storing the output data, user can do the followings: 

1. Always keep an original file of the Microsoft '' Excel code. 

2. Save the file with a new filename. 

3. From the menu bar, Select File->Save As and name the file, e.g., offdesign-051704.xls. 

4. Click Save button to confirm saving. 

5.4.2 Troubleshooting . — This program has been debugged several times. Therefore whenever you 
encounter computation errors due to either zero division or square root of a negative quantity, you will be 
notified through a pop up window indicating where the computation problem is. If you have any 
comments or bug problems you encounter in the program, you can report them to us for further 
improvement. Detail about the contact number can be found in the ‘CoverPage’ sheet of the program. 

In multiple -point calculation, users may always see some of these pop-up windows as shown in 
figure 5.5(a) to (d). 

There are three possibilities of having this error or warning messages: 

(i) During the multiple-point computation, the design variable may be out of the range. It will 
happen, especially when the design variable such as flight Mach number is out of the range. For instance, 
the thrust of an engine is approaching zero as the flight Mach number is 5.1. If the upper limit of flight 
Mach number is set to be greater than 5.1 (e.g., 6 ), the code will not be able to proceed at Mach number 
higher than 5. 1 and thus give the warning messages shown in figure 5.5(c). 

(ii) The iteration is not converging for some engine configurations at the specified flight Mach 
number and altitude. 

(iii) The input values are incorrect. For instance, the main burner exit temperature ( T t 4 ) may be too 
small for an engine to produce enough thrust. Please double check all input values and try again. 



(a) Example of warning messages 



(c) Equation where the error occurs 

Figure 5.5. — Pop-up windows 



(c) Example of warning messages 



(a) Current values of flight condition 
showing warning messages 
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5.5 Summary 


In summary, user needs to do the followings to run the program: 

Go to the “Input” sheet. 

Specify the Unit System (SI or English?) 

Specify the ITB Switch. (Is the computation for engine with ITB-ON or for engine with 
ITB-OFF?) 

Specify the Fuel Burn Model. ( CSH or MSH model?) 

Import Reference Conditions data from on-design code. 

Enter all the input parameters indicated in ‘dark green’ cells (do not modify or change the value 
indicated in ‘hotpink’ cells). 

Specify Perform Computation in. User will be directed to a new control sheet depending on the 
selection. 

If you select “ Full Throttle Performance (FTP): Multiple-point do the followings in 
‘ControlFTP’ sheet: 

Specify a range for flight Mach number (Mo). 

Specify number of altitude (new ‘dark green’ cells will be created). 

Specify all input parameters indicated in ‘dark green’. 

Click Calculate button to compute and plot the results. 

Repeat the above steps for different input parameters. 

Go to the ‘Input’ sheet if you want to change some engine input parameters. 

If you select “ Full Throttle Performance (FTP): Single-pointf do the followings in ‘SinglePt FTP’ 
sheet: 

Enter a value in flight Mach number, (Mo) cell. 

Enter a value in Altitude (h) cell. 

Click Calculate button to compute and display the results. 

If you select “ Partial Throttle Performance (PTP)f do the followings: 

a. Specify a value for altitude (h). 

b. Enter a value in Min % thrust cell. 

c. Specify number of flight Mach number (new ‘dark green’ cells will be created). 

d. Specify all input parameters indicated in ‘dark green’. 

e. Click Calculate button to compute and plot the result. 

f. Repeat the above steps for different input parameters. 

g. Go to the ‘Input’ sheet if you want to change some engine input parameters. 
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Appendix A — U.S. Standard Atmosphere 1976 


A.l Introduction 

The US Standard Atmosphere 1976 is an average, piece-wise continuous function for density, 
pressure and temperature profiles of the earth’s standard atmosphere. The standard atmosphere is 
mathematically defined in seven layers from sea level to 84.852km. It is a function of geo-potential 
altitude, which is represented by: 


^ ' R earth 
h + R ea rth 


(A.l) 


where h - geometric altitude [km] 

Rearth ~ radius of the earth [= 6356.66km] 

The table below summarizes the values of several parameters at each of the defined levels. 


Table At. — Values of several parameters at each defined levels of the earth’s atmosphere. 


Layer 

(0 

Name 

Lower 
Geopotential 
Altitude (z) 
[km] 

Upper 

Geopotential 

Altitude 

[km] 

Lapse rate 

(LR st d,i) 

[K/krnj 

Standard 
Pressure, (P st d,i) 
[Pa] 

Standard 
Temperature, 
( T s td,i ) 

[K] 

1 

Troposphere 

0 

11 

-6.5 

101325 

288.15 

2 

Stratosphere 

11 

20 

0.0 

22632.06 

216.65 

3 

- 

20 

32 

+1.0 

5474.89 

216.65 

4 

- 

32 

47 

+2.8 

868.02 

228.65 

5 

- 

47 

51 

0.0 

110.91 

270.65 

6 

Mesosphere 

51 

71 

-2.8 

66.94 

270.65 

7 

- 

71 

84.852 

-2.0 

3.96 

214.65 


A.2 Temperature 

The variation of temperature, T(z ), with geo-potential altitude is expressed by a continuous and 
piecewise linear relation. If Lapse rate is non-zero (LR st dj ^ 0), 

T ( z ) = T std,i + LRstdd ■ ( z - z std,i ) (A.2) 


If Lapse rate is zero ( LR stc ij = 0), 


R( z ) = R s td,i 


(A-3) 


A.3 Pressure 

The pressure, P(z), is also a continuous piecewise function of geo-potential altitude. If Lapse rate is 
non-zero (LR stdi V 0), 
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If Lapse rate is zero ( LR stc a = 0), 


f 

P(z) = P st d,i ■ ex P 

V 


-1000 ■ g 0 -W 0 • (z ■ 


R -T 


std,i 


where 

go - acceleration of gravity (= 9.80665 m/s 2 ) 

R* - gas constant for air (= 83 14.32 J/kmol-K) 

Wq - Molecular weight of air (= 28.9644 kg/kmol) 
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Appendix B — Combustion Fuel Burn Models 

B.l Types of Fuel Burn Models 

Two models (p. 116 of ref. 2) are being used to calculate the fuel burned in the main burner and the 
ITB: 

(1) Constant Specific Heat ( CSH) model — The air and combustion gases at inlet and exit of each 
component are modeled as “calorically perfect gases” with constant specific heats. The values of the 
specific heats are different at inlet and exit of two combustors (main burner and ITB) 

(2) Modified Specific Heat ( MSH) model — All engines properties are calculated using CSH model 

except the fuel used. In this model, the inlet total temperature (T t j and of combustors are calculated 

using CSH model while the exit total temperatures (T l4 and T t4 f) are directly obtained from the user 
inputs. Nevertheless, the total enthalpies (h t4 , h t4 , h t4 4 , and h t4 f) will be calculated directly from the 
Variable Specific Heat (VSH) model (p. 1 16 of ref. 2). Provided hpR is given, the amount of fuel burned 
is then calculated using eq. (B4) and (B7), respectively. Therefore, the improvement on enthalpy 
calculation gives better estimates of fuel used. 

The following ENTHALPY subroutine is based on pages 105 and 106 of reference 3 and used in the 
MSH model to calculate the enthalpy (h) at engine stations t3, t4, t4.4, and t4.5. 


Subroutine Enthalpy 

ENTHALPY (T,f h) 

Inputs: T and / 

Output: h 


hnir 


h ref + a 0 t + A t 2 + A_ r 3 + ^-r 4 + A. t 5 +^t 6 + A. t 1 + ^-r 8 

2 3 4 5 6 7 8 


(B.l) 


h 


prod ~ 


Kef + A 0 T + A T 2 + A r 3 + A J4 


+ ^±T 5 +^T 6 +^T 7 + 
5 6 7 



(B-2) 


hair fhprod 


1 + / 


(B-3) 


where the constants Aq, etc. are given in table B 1 (p. 1 06 of ref. 3). 
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Table Bl. — Constants for air and combustion products used in subroutine ENTHALPY(0, T,f h, c p , y). 


Constant 

Air alone 

Combustion products of air and (CIDn fuels 

A 0 

2.5020051 x 10“' 

7.3816638 x 10 _1 

Ai 

-5.1536879 x 10“ 5 

1.2258630 x 10“ 5 

a 2 

6.5519486 x 10' 8 

-1.3771902 x 10' 8 

A 3 

-6.7178376 x 10' 12 

9.9686793 x 10“ 12 

a 4 

-1.5128259 x 10“ 14 

-4.2051104 x 10“ 14 

As 

7.6215767 x 10“ 18 

1 8 

1.0212913 x 10 

A 6 

-1.4526770 x 10“ 21 

-1.3335668 x 10“ 21 

Ay 

1.0115540 x 10“ 25 

7.2678710 x 10“ 25 

href 

-1.7558886 Btu/lbm 

30.58153 Btu/lbm 


The equations needed to calculate the fuel-air ratios of the main burner (//,) and the ITB (//,/,) are listed 
below for both the CSH and MSH models. 

Main burner 

If CSH model then (eq. (27) of ref. 5) 


fb = 


XX-b XyXcLXcH 
r \bhPR/ho~XX-b 


(B-4) 


Else 


T t3 - To X,X C LX C H 


(B.5) 


ENTHALPY (0,T t3 ,0,h t3 ,0,0) 

Set initial value of fuel/air ratio at station 4 


3 ENTHALPY (0,r t4 ,f 4h k t4 ,0,0) 


y _ ht4~h t3 

r\bhpR-h t 4 

If \ fb - f/\j | > 0.0001, then j\j = fo and go to 3; else continue. 
End if 

Interstage Turbine Burner 

If CSH model then (eq. (30) of ref. 5) 


fub={\ + fb) 


t X-itb ~ tX-b^tH 
r\i tbhpR/ho -t x-itb 


(B-6) 


(B-7) 
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Else 


Tt4A~Tt4ltH 

ENTHALPY {0,T t4A ,f b ,h t4A ,0,0) 

Set initial value of fuel/air ratio at station 4.4 =f 4 .s\ 

4 ENTHALPY {0,T t4 . 5 ,f 4 . 5i ,h t4 . 5 ,0,0) 

fitb = (1 + fb) hl4 - 5 ~ ht4A (B.8) 

r \it b hpR-h t4 .5 

If | f b + f itb - / 4 5; j > 0.0001, then / 4 5j = f b + f ub and go to 4; else continue. 

End if 

lf fitb = 0 then Yi* = Y r c pitb =C pt , and R itb - R, 
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Appendix C — Control Limits 


C.l Introduction 

Normally, the designers of an aircraft gas turbine engine need to consider the engine control system 
so that the newly-designed or improved gas turbine engine operates well within its flight envelope. For 
example, the throttle setting needs to be controlled so that the burner exit temperature will never exceed 
the turbine blade maximum temperature. Another example is that fan or compressor surges or stalls 
should not occur during the operation. 

Therefore, engine’s control system plays a significant role in this engine performance cycle analysis 
because of the several advantages: (1) maintain stable airflow, internal pressures and temperatures, and 
rotor speeds within safe operating limits; (2) avoid significant speed, pressure or temperature variation, 
and stalls or surges (ref. 2). In addition, it also demonstrates how an aircraft engine behaves away from 
their reference or on-design conditions. 

The engine’s control system in this study is specifically designed for a separate-flow, dual-spool, 
uncooled turbofan engine with ITB. Moreover, the development of the equations is based on the 
foundational materials found in references 1 and 2. 

C.2 Variables of Engine Control System 

C.2.1 Dimensionless Freestream Total Temperature . — Due to its significant influence on overall 
behavior of the turbofan engine, a dimensionless form of the freestream total temperature is defined, i.e., 
the ratio of the freestream total temperature to the sea-level static ( SLS) temperature of the standard 
atmosphere, as shown below (ref. 2): 


On = 


T t o 


1 + 


Jc -1 


M- 


T 

1 r 


ref 


T, 


- = T, 


ref 


3 _ 

T re f 


(C.l) 


(Note:T,. e/ =T std ) 

This property is called theta 0 (Go) and is used to combine two parameters, i.e., altitude (To) and flight 
condition (Mo) into a single parameter. Therefore, every point in the flight envelope has a specific value 

of 00- 

In addition, 9o has the following properties: 

(1) Go = 1 at sea level static conditions; 

(2) 0o can be greater or less than 1.0; 

(3) 9o depends only on Mach number above between 36.5 and 65 kft (11.12 and 19.81 km). 

C.2.2 Throttle Ratio . — The other useful property of an engine control system is the throttle ratio 
(77?). It is defined as the ratio of the maximum T t 4 value to the T t 4 value at SLS conditions. It can be 
expressed as: 

77? = r?4max (C.2) 

TfASLS 

When the reference altitude is at SLS ( h=0 ), then T^sls — T t 4 R ■ Otherwise, TR has to be manually 
input by user. 
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Figure Cl. — Compressor Ratio as a Function of 0o and 7 ' /4 of a Single-Spool Turbojet 
Engine for a Compressor with a Reference Point of n c - n c max = 20 , 

T,A = T t 4 max = 3300 °R , and 0 O = 1.1 (ref. 2). 


C.2.3 Theta Break . — Normally, the change in shape of the performance curve occurs at simultaneous 
maximum of n c and 7' /4 . In figure 3. 1 , it is clear that the control logic must switch from limiting n c to 
limiting t i4 at this point of maximum n c and T t4 . This unique point is known as the theta break 
(00 break)- 


C.3 Types of Control l imits 

In this study, only three types of engine control limits will be considered, i.e., turbine blade maximum 
temperature, compressor limitations and engine speed. 

C.3.1 Material Limitation . — It is known that an increase in permissible turbine inlet temperature will 
lead to a reduced TSFC and increased specific thrust (ref. 4, p. 357). Although the current advanced 
technology in material allows higher inlet turbine temperature, it is still essential to limit T tA or 7 ) 4 5 for 
safety purpose. 

C.3. 2 Compressor Pressure Ratio. — Three limits, namely, 7 t cmax , P, 4 max and 7) 3max , are selected to 
stabilize the internal behavior of the compressors and avoid compressor surges or stalls. 

C.3.3 Engine Speed . — As shown in section 5.2.7 of reference 2, the change in total enthalpy across a 
compressor or a fan is proportional to the rotational speed N squared for a calorically perfect gas. 

C.3. 3.1 Low-pressure Spook For a low-pressure compressor ( LPC ), it can be expressed as follow: 

h t 2. 5 -h t 2=K\N 2 ip (C.3) 
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Using referencing, it can be rewritten as: 


Nlp I h t2 '5 - h^~ _ I h t2 {x c L -\ ) 

N LPR v htl.SR - h t 2R v ht2R (?cLR ~ l) 

I h t o(x c L-l) „ I T, q{x cL - l) 

V KoR^cLR -l) V TtOR^cLR _ l) 

9o(x c z -!)" 

V 00 r{i C LR -l) 


(C.4) 


Since the LPC and the fan are on the same shaft, the enthalpy rise across the LPC will be proportional to 
the enthalpy rise across the fan during normal operation, (ref. 2) 


X cL "I 
X f -1 


f P\ 
*cL “I 


C / 


-1 


(C.5) 


Then, equation (C.4) can be rewritten as 


N LP Go (V 1) 

N LPR \6or( x JR-1) 


C.3.3.2 High-pressure Spool. For high-pressure compressor ( HPC ), it is in the similar form: 


h t2 


— h,2.5 = K^N 


2 

HP 


Using referencing, it can be rewritten as: 

N HP I hfi ~ ht2~ _ I ll t2.5^cH ~l) 

N I IPR V ht3 R ~ h t 2.5R v ht2.5R (t cHR ~ l) 

I h t QXcL{^cH ~l) ~ I T t QX c L{^ c H ~ l) 

V htORtcLR (x cHR - 1) V T tORtcLR (x C HR ~ 0 

N HP 0Q x cl{ t cH ~l) 

N IIPR V 00 RlcLR iScHR - l) 


(C.7) 


(C.8) 


(C.9) 


C.4 Relationship between T t 4. 5 and n c 

C.4.1 High-pressure compressor . — From a power balance between HPT and HPC, one will obtain 
the performance equation for the total temperature ratio of the HPT. 


x cH - l + r l mH ( 1 + /&) ~ b — 

I,.X C £ 


(C.10) 
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T c tf=l + C 0 ^— (C.15) 

OO x cL 

Substituting (C. 1 5) into (C. 1 1 ) gives 

Orff-jl + C,5 £ J-)}l , A-') (C.16) 

l 00 *cL j 

where 

C, = C 0 r\ cH (ii) 

Compared to equation (D.3) of reference2 (pp525), k c h is not only a function of throttle setting ( T t4 ) 
and flight conditions (Oq), but also a function of x c j : 
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C.4.2 Low-pressure compressor 


From a power balance between LPC, fan and LPT, one will obtain the performance equation for the 
total temperature ratio of the LPT. 


a ( x / “ 0 + ( IcL - l) = 11 mL (l + fb+ fitb (l - T tL ) 

T r 

Substituting equation (C.5) into (C.18), it yields: 

v = 1 + (yn - ih„t lhzML ( ' ~ - y + j - 1 

The LPC total temperature ratio follows from equation (C.5): 


X cL 




Using LPC efficiency, the LPC total pressure ratio is given by: 

^cL = {l + ri c t(x £ i-l)} Yc ^ Yc ^ 

From equation (6) of reference 4, 

C pitbT t4.5 

X X—itb ~ 

Substituting equation (22) into equation (19) yields 

CpitbTfA.S 


"/ + - ')iw ■ 


Cpc-To 


C pc T o 


( T 0 } 

[ (l _ l tL X 1 + fb + fitb )| 

JrefQo , 

^cLR -l + a(lyR -l) 


(C.18) 


(C.19) 


(C.20) 


(C.21) 


(C.22) 


(C.23) 


Similar to the assumption in HPT, the flow is choked at station 4.5. Therefore, the value of x t i will 
remain constant. In addition, //, and fpp are negligible when compared to unity. 

Unfortunately, a is not a constant but a function of (n c i , n c fj Let 


( , 5 


(l-X^Xl + fb+fitb) 

L pc ytref J 


(iii) 
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Then, 


Ty — 1 + C 2 


T, 4.5 f 1 

9 0 \t cLR - l + a(x f R -l) 


Substituting (C.24) into (C.20) yields 


t c l -I + C 3 


T t 4.5 f 1 

9 0 [t c z,k - l + a(x^ - l) 


(C.24) 


(C.25) 


where 


f-3 — ^ ]mL 


r 

^ rm 


pitb 


c 


pc 


l re f 


(l ~1tL )(l + fb + fitb \lcLR ~ l) 


Substituting (C.25) into (C.21) yields 


ft cL 


l+a 


T t 4.5 


ftcLR -l + a(tjK -l) 


YcAY e _1) 


(iv) 


(C.26) 


where 


Q -TciQ 


(V) 


Again, when compared to equation (D.3) of reference 2 (pp525), n c p is not only a function of throttle 
setting (T t f) and flight conditions (9 q), but also a function of a : 


f 

ftcL = J 

V 


Tt 4.5 

9o 


\ 

,a 

/ 


(C.27) 


C.5 Limiting procedures 

C.5.1 Reference conditions . — The values of T t4 max and Tt4.5 max (control limit) can only be smaller 
or equal to the values of T t4R and T t45R at the reference ( on-design ) condition. Furthermore, when the 
reference condition is at a reduced Tq R due to the altitude, the reference T t4R and T t45R have to be 
reduced too. This may be followed by the reduction in the values of T t4 max and T t 4.5 max too. 

%r=t, r ^ (C.28) 

‘ ref 

If 0o« > TR then 

Tt 4 max - T t4R 
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else 


T t4 .5 max — 7)4.5/? 


(C.29a) 


T, 4 max — Tt4R 
7)4.5 max — 7)4.5/? 


' %R _ A 

v TR 


OR 

V TR 


f 0 n A 


(C.29b) 


C. 5.2 Limiting 7) 4 and 7/4.5 . — As shown in appendix D of reference 2 (p. 525), as 9o is decreased to 
the left of the theta break, or 9o break 7)4 must be reduced in order to protect the compressor from poor 

T 

operation. By doing so, the term — — must remains constant. 

9o 

A simple limiting procedure as shown below is included into the code in order to limit 7)4 values 
equal to or less than 7 ) 4 max . 


9o = t T 


Jo_ 

Tref 


If 9 0 > TR then 


T t 4 


- 7/4 1 


else 


7)4 


_ T I 0 Q 

— 1 t4 max | 


(C.31) 


(C.32) 


(C.33) 


When ITB is considered, a separate limiting procedure is required to limit 7) 4 5 values equal to or less 
than 7)4 5 max . 

If9 0 >r7? then 


7)4.5 = 7 ) 4.5 max (C.34) 

else 


7)4.5 =7) 4 .5max(jy (C-35) 

C5.i Limiting n f and ti c u . — For a two-spool turbofan engine with ITB, there are two engine fuel 
flows. The fuel flow to the main burner directly influences the speed of the HP spool and the HP 
compressor pressure ratio. The fuel flow to the ITB directly influences the speed of the LP spool and the 
fan and LP compressor pressure ratios. 
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In Full Throttle Performance, when HP compressor pressure ratio is greater than its given maximum, 
the actual engine control will reduce the required fuel flow to the main burner in real time until the given 
limit is reached. In this study, the code does not calculate the reduction in fuel flow directly. Instead, the 
code will reduce T t 4 (main burner) or T t4 5 ( ITB ) to a lower value. 

Therefore, two control limits have been selected in Full Throttle Performance, i.e., fan pressure ratio, 
n / max (LP spool) and HP compressor pressure ratio n c Hmax (HP spool). Based on this, two throttle 
control loops have been added to the code as shown in table Cl. 

When the fan pressure ratio ( FPR ) is greater than its limit, n / max (case 1), the control loop will 
reduce T t 4 5 (reduce fuel flow to the ITB). When the HP compressor pressure ratio ( HPC PR) is greater 
than its limit, k c h max (case 2), the control loop will reduce T t4 (reduce fuel flow to the main burner). 

If ITB is on, the control loops for FPR and HPC PR are sequential, i.e., FPR limit control goes first, 
followed by HPC PR limit control. 

If the ITB is off and the FPR is still greater than its limit, n / max (case 3), the control loop will reduce 
T t 4 (reduce fuel flow to the main burner). This reduction in T t4 will reduce both FPR and HPC PR. 
Similarly, this reduction in T l4 applies to a case where the FPR is still too high even when 7 ) 4 5 is 
reduced until ITB is off (case 1 shifted to 3). 


Table Cl. — The logic sequence of reducing temperature (reducing fuel flow to the burners) 


Parameter 

‘reduce’ 

Reduce T t 4 . 5 when 

Reduce T t 4 when 

Case 

max 

71 cH - n cH max 

n f - ft/ max 

k cH - KcHmax 

ITB on 

1 

- 

- 

2 

ITB off 

- 

- 

3 

4 
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Figure C2. — Flow-Chart of Control Limits for Fan Pressure Ratio 
and High-Pressure Compressor Pressure Ratio. 
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Appendix D — Flowchart of Iterative Solution Scheme 



Figure Dl. — Flow-Chart of the Iterative Solution Scheme 
for Full-Throttle Performance. 
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Figure D2. — Flow-Chart of the Iterative Solution Scheme 
for Partial-Throttle Performance. 
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Figure D3. — Flow-Chart of the Iterative Solution Scheme for 
the Subroutine itbperfl (Central ITB Computation) 
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Appendix E — Reference Engine Data From On-Design Code 


E.l Introduction 

In general, on-design (parametric) cycle analysis differs significantly from off-design (performance) 
cycle analysis. Thus, two separate codes using Micro so ft" Excel and VBA are written for each cycle 
analysis. 

Once the most promising engine is selected by user (i.e., all engine components are fixed, e.g., CPR, 
BPR, altitude, FPR etc.), it is then run using on-design code to get the engine performance at the selected 
operating conditions. 

These output data represents the selected engine performance at a design-point. In order to predict the 
engine performance at off-design conditions, these outputs are then used as the reference conditions in 
off-design code. This step is important because different reference conditions (i.e., different engines) will 
give different performance over the whole flight spectrum. 

For convenience, authors have modified the on-design code to include 'Input/Output files' features in 
each code. By doing so, the output from on-design code can be saved as the reference engine data file, 
which is required in off-design calculation. As a result, users do not need to manually input the reference 
conditions, which is a time-consuming process. 


E.2 Instructions to get reference engine data 

Follow the step 1 to 12 to generate a reference engine data file. 

1. Generate the sample reference engine data as shown in table El using on-design code (ref. 4). 

2. Open off-design code. 

3. Select ‘Input’ sheet as shown in figure 5.1. 

4. Click on ‘ Import data from ITB ondesign code ’ button. 

5. A dialog box similar to figure El will appear and request for specifying the path and filename. 

6. The default path and filename is “c Preference, txt”. Users may need to enter a new path if 
necessary, for example, “ c:\ITB\reference.txf\ 

7. An input dialog box as shown in figure E2 appears when one of these situations happens: 

(7a) The specified directory (default drive is 7 i:V) does not exist. 

(7b) The file (e.g., “ refelrence.txt ” in figure E2) does not exist. 

Check for existing drive path or/and specify another filename. 

8. A dialog box as shown in figure E3 appears when the unit system in reference data file is 
different from the current unit system in off-design code. Change the unit system in ‘Unit 
System’ pull-down list. Then, start from step 2. 

9. An input dialog box as shown in figure E4 appears and request for entering the value of mass 
flow rate. The default value is 100 kg/s. 

10. Enter value and click 'Ok' button to proceed. 

1 1 . All reference conditions data and engine characteristics of a specific design-point engine are 
successfully imported. The specified pathname will be shown in ‘Reference-point engine path:” 

12. It is ready to perform off-design calculation. 

13. Repeat step 1 to 10 for different engines. 
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Table El. — Input of Reference Data for on-design code 


Reference Conditions 

Input 

Mach number ( Mor ) 

0 

Altitude (h, kft) 

0 

Main burner exit total temperature (T/ 4 r, °R) 

3600 

ITB exit temperature (A/ j^, °R) 

3275 

LPC pressure ratio (k c lr) 

2.31 

Compressor pressure ratio (n c j^) 

26.8 

Fan pressure ratio (ji/k) 

2.31 

Fan bypass ratio (a r) 

1.91 

Mass flow rate ( m 0R , lbm/s) 

1097 


ITB cycle analysis 


Enter the path and filename (that is not used by 
another program) 



OK 


Cancel 


Ih: \reference.txt 


•igure El. — Dialog Box for Specifying the Path and Filename 
where the Reference Conditions Data Will Be Stored. 



Figure E2. — Dialog Box when there is An Error Reading Reference 
Conditions Data from the Specified Path and Filename. 


ITB performance cycle analysis 


Unit system inconsistent! Please change unit and try again. 


OK 


Figure E3. — Dialog Box when there is Unit System Inconsistency. 
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Figure E4. — Dialog Box for Inputting 
Value of Mass Flow Rate. 
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Appendix F — Visual Basic Codes 

The Visual Basic codes can be easily seen on the spreadsheet under “ToolsYMacroYVisual Basic 
Editor”. 


Sub itbperflQ 
Dim etar 
Style = vbYes 

Title = "ITB Performance Cycle Analysis" 

Stylel =vbYesNo 

t 

If plottype = 1 Then 

msgO = "Occurs at " & "MO = " & MO & ", h =" & alt & 

Else 

msgO = "Occurs at " & "MO = " & MO & ", Tt4 =" & Tt4 
End If 
bigit = vbNo 

I 

'Set initial values for some engine components 
mdotO = mdotOR 
fb = f_bR 
fitb = fitbR 
pitH = pitHR 
pitL = pitLR 
picH = picHR 
picL = picLR 
pi_f = pi_fR 
tautH = tautHR 
tautL = tautLR 
taucH = taucHR 
taucL = taucLR 
tauf = taufR 

f 

'Restore values for some parameters if lTBswitch is internally ON again in next MO iteration 
If ITBswitch = "ON" Then 
tau itbR = Tt4p5R / (Tt4R * tau tHR) 'eqn(5r) 

'cp itb = Sheets("Input").Cells(7, 14) * conversion 'removed 3/8/05 
'gama_itb = Sheets("Input").Cells(10, 14) ' 

Else 

tau_itbR =1# 

End If 

t 

'Preliminary computations 
tau_r = 1# + cnst5 * MO * MO 'eqn(l) 
pi_r = tau r A cnst8 'eqn(2) 

If MO <= 1# Then 

eta_r=l# 'eqn(3) 

End If 

If MO > 1# And MO <= 5 Then 
eta_r = 1# - 0.075 * (MO - 1#) A 1.35 
End If 

pi_d = pi dmax * eta r 'eqn(4) 

'R itb = (gama_itb - 1#) * cp_itb / gama_itb 'eqn(7) 
aO = (gama_c * R_c * gc * TO) A (0.5) 'eqn(8) 
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MFPMO = MO * cnst4 * taur A cnst6 'eqn( 1 1) 

TtO = TO * tau r 'eqn( 10) 

PtO = PO * tau r A cnst8 'eqn(9) 

itblink = 1 
10 

End Sub 
Sub itbperCQ 
Dim pi_tHN, taufN 

Dim MFPM8, MFPM9, MFPM18, MFPM19 

Dim ebrl, ebr2, flpc 

Dim h_t3, f_4i, Tt4p4, h_t4p4, f_4p5i 

Title = "ITB Performance Cycle Analysis" 

Style = vbYes 
Stylel =vbYesNo 

t 

h_0 = cp_c * TO 
If h_0 = 0# Then 

Msgl 1 = " Division by zero in 'h_0 equation'. " 
response = MsgBox(Msgl 1, Style, Title) 

Msgl 1 = " Check if' TO = O'" 
response = MsgBox(Msgl 1, Style, Title) 

Msgll = " Current values => cpc =" & cp_c & ", T0=" & TO 
response = MsgBox(Msgl 1, Style, Title) 
response = MsgBox(msgO, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Stylel, Title) 
itblink = 0 
GoTo 10 
End If 

limit = 0 
loopif = 0 
loopich = 0 
loop 1 = 0 
reduce = 1 
pi cHo = pi cHR 
taucHo = taucHR 
GoTo 19 
'liew = 1 

'Check intermediate values for debugging 
Sheets("test").Cells(l, 1) = "MO" 

Sheets("test").Cells( 1 , 4) = "loop2" 

Sheets("test").Cells(l, 2) = "loopif' 

Sheets("test").Cells( 1 , 3) = "loopicH" 

Sheets("test").Cells( 1 , 5) = "Tt4" 

Sheets("test").Cells( 1 , 6) = "Tt4p5" 

Sheets("test").Cells( 1 , 7) = "tau_fN" 

Sheets("test").Cells(l, 8) = "tau_cH" 

Sheets("test").Cells( 1 , 9) = "pi_f' 

Sheets("test").Cells( 1 , 10) = "pi_cH" 

Sheets("test").Cells(l, 11) = "overpi_f' 

Sheets("test").Cells(l, 12) = "overpi_cH" 

Sheets("test").Cells( 1 , 13) = "ITB" 

Sheets("test").Cells(l, 14) = "reduce" 

19 

If ITBswitch = "OFF" Then reduce = 3 
12 
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GoTo 29 

'Check intermediate values for debugging 
If loopich > 0 Then 
Sheets("test").Cells(liew +1,1) = MO 
Sheets("test").Cells(liew +1,4) = loop2 
Sheets("test").Cells(liew +1,2) = loopif 
Sheets("test").Cells(liew +1,3) = loopich 
Sheets("test").Cells(liew +1,5) = Tt4 
Sheets("test").Cells(liew +1,6) = Tt4p5 
Sheets("test").Cells(liew +1,7)= tau_f 
Sheets("test").Cells(liew +1,8)= tau_cH 
Sheets("test").Cells(liew +1,9)= pi_f 
Sheets("test").Cells(liew + 1, 10) = pi_cH 
Sheets("test").Cells(liew + 1, 11) =overpi_f 
Sheets("test").Cells(liew + 1, 12) =overpi_cH 
Sheets("test").Cells(liew + 1, 13) =ITBswitch 
Sheets("test").Cells(liew + 1, 14) = reduce 
liew = hew + 1 
End If 
29 

'Recalculate cp_t, gama t and R_t after either (i) controlimit or (ii) Tt4 cut 
Call enthalpy(l, Tt4, f_b, 0, cp_t, gama_t) '3/8/05 
cnst9 = (gama_t - 1#) / gama_t '3/8/05 
R_t = cnst9 * cp_t '3/8/05 

tau_lb = cp_t * Tt4 / h_0 'eqn(15) 

t 

'Component: INTERSTAGE-TURBINE BURNER 
1 

If loop 1 > 20 Then 
loopl = 21 
itblink = 0 
GoTo 10 
End If 

11 

'Recalculate cp itb, gama itb and R itb after either (i) controlimit or (ii) Tt4 cut 
Call enthalpy) 1, Tt4p5, f_itb, 0, cp_itb, gama_itb) '3/8/05 
R_itb = (gama_itb - 1#) / gama_itb * cp_itb '3/8/05 

loop2 = 0 

If ITBswitch = "ON" Then itbinit = 1 'added at 6/1 1/04 
tau_litb = cp_itb * Tt4p5 / h_0 'eqn(16b) 'added at 6/1 1/04 

'Check if ITB is still operating 

If tau_litb <= (tau_lb * tau_tH) Then itbinit = 0 'added at 6/1 1/04 
If ITBswitch = "ON" And itbinit = 0 And reduce = 1 Then 'added at 6/1 1/04 
'If ITBswitch = "ON" And Tt4p5 / (Tt4 * tau tH) <= 1 And reduce = 1 Then 
ITBswitch = "OFF" 
loopif = 0 
reduce = 3 
End If 

'Make changes to some parameters while 1TB is turned off 
If ITBswitch = "ON" And reduce = 1 Or reduce = 2 Then 
tau_itb = Tt4p5 / (Tt4 * tau_tFl) 'eqn(16a) 

' tau_litb = cp_itb * Tt4p5 / h_0 'eqn(16b)'removed at 6/1 1/04 
Else 

tau_itb=l# 'eqn(16c) 
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tau_itbR = 1# 

tau_litb = cp_t * Tt4 * tau_tH / h_0 'eqn(16d), ignore tau_itb=l 

Tt4p5 = Tt4 * tau_tH 'eqn(16e) 

Call enthalpy(l, Tt4p5, 0, 0, cp_itb, gama_itb) '3/8/05 
R_itb = (gama_itb - 1#) / gama_itb * cp_itb '3/8/05 
End If 

A4p5oR = (tauitb / tau_itbR) A n_a 'eqn(17) 

A8oR = (tau itb / tauitbR) A m_a 'eqn(18) 

1 

'Component: HIGH-PRESSURE TURBINE 
tau_tH = 1# - eta_tH * (1# - pi_tH A cnst9) 'eqn(19) 

pi_tHN = (tau_tH * tau_itb / (tau_tHR * tau_itbR)) A (0.5) / A4p5oR * pi_tHR 'eqn(20) 
If Abs(pi_tHN - pi tH) > 0.00001 Then 
pitH = pi_tHN 
loopl = loopl + 1 
GoTo 1 
End If 

2 

If loop2 > 50 Then 

Msg 17 = "Eqn(38) => 'tau_fN' does not converge!" 
response = MsgBox(Msgl7, Style, Title) 

Msgl7 = "Reduce or increase 'pi_f_max' value and try again." 
response = MsgBox(Msgl7, Style, Title) 
response = MsgBox(msg0, Style, Title) 

If singlept = False Then jump = MsgBox(MsgJump, Stylel, Title) 
itblink = 0 
GoTo 10 
End If 

t 

'Component: HIGH-PRESSURE COMPRESSOR 
tau_cH = 1# + eta_mH * (1# + f_b) * tau_lb * (1# - tau_tH) / (taur * tau_cL) 'eqn(21) 
pi_cH = (1# + eta_cH * (tau_cH - 1#)) A cnst8 'eqn(22) 
err_picH = Abs(pi_cH - picHmax) 

If limit = 4 And err_picH > 0.001 Then 
DpicH = pi cHo - pi_cH 

If Dpi cH = 0 Then Dpi cH = 0.00001 'avoid division by zero 
DTt4 = Tt4o - Tt4 
pi_cHo = pi_cH 
Tt4o = Tt4 

DifTt4 = (pi cH - pi cHmax) * DTt4 / Dpi cH 

If Abs(DifTt4) > 100# Then DifTt4 = DifTt4 / Abs(DifTt4) * 100# 

Tt4 = Tt4 - DifTt4 
limit = 4 
GoTo 12 
End If 

'reduce = 2 Or reduce = 4 Then 
'skip control limit when PTP 

f 

'Component: MAIN BURNER 
temp_fb = eta_b * h_PR / h_0 - tau_lb 
If temp fb = 0# Then 

Msgl2 = " Division by zero => Refer to eqn(21a)." 
response = MsgBox(Msgl2, Style, Title) 

Msg 12 = "Eqn(21a) => Check if ' eta b * h_PR / h_0 = tau lb '" 
response = MsgBox(Msgl2, Style, Title) 
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Msg 12 = "Current values => eta_b=" & eta b & ", h_PR=" & h_PR & ", h_0=" & h_0 & ", tau_lb=" & 

taulb 

response = MsgBox(Msgl2, Style, Title) 
response = MsgBox(msgO, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Stylel, Title) 
itblink = 0 
GoTo 10 
End If 

'Fuel model 1: CSH 
If fuelmodelname = "CSFI" Then 
f_b = (tau_lb - taur * tau_cL * tau_cFl) / temp_fb 'eqn(23a) 

Else 

'Fuel model 2: MSH 
loop3 = 0 

Tt3 = TO * tau r * tau_cL * tau_cFI 'eqn(23b) 

'Call enthalpy(Tt3, 0, h_t3)'3/8/05 
Call enthalpy(0, Tt3, 0, h_t3, 0, 0) 

f_4i = (tau_lb - tau r * tau_cL * tau_cFl) / temp_fb 'eqn(23c)'initial value for f_4i 
3 

If loop3 > 20 Then 
itblink = 0 
GoTo 10 
End If 

'Call enthalpy(Tt4, f_4i, h_t4)'3/8/05 
Call enthalpy(0, Tt4, f_4i, h_t4, 0, 0) 
temp_fb = eta_b * h_PR - h_t4 
If temp fb = 0# Then 

Msgl2 = "Division by zero => Refer to eqn(21d)." 
response = MsgBox(Msgl2, Style, Title) 

Msg 12 = " Check if' etab * h_PR = h_T4 '" 
response = MsgBox(Msgl2, Style, Title) 

Msg 12 = "Current values => eta_b=" & eta b & ", h_PR=" & h_PR & ", h_T4=" & h_t4 
response = MsgBox(Msgl2, Style, Title) 
response = MsgBox(msg0, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Stylel, Title) 
itblink = 0 
GoTo 10 
End If 

f_b = (h_t4 - h_t3) / temp_fb 'eqn(23d) 

If Abs(f_b - f_4i) > 0.0001 Then 
f_4i = f_b 
loop3 = loop3 + 1 
GoTo 3 
End If 
End If 

If ITBswitch = "OFF" Then 'skip fitb computation 
Call enthalpy(l, Tt4, f_b, 0, cp_t, gama_t) '3/8/05 
cnst9 = (gama_t - 1#) / gama_t '3/8/05 
R_t = cnst9 * cp_t '3/8/05 
f_itb = 0# 

GoTo 5 
End If 

f 

'Component: INTERSTAGE-TURBINE BURNER 
temp_fitb = eta_itb * h_PR / h_0 - tau_litb 
If tempfitb = 0# Then 
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Msgl3 = "Division by zero => Refer to eqn(22a). Check if 'eta itb * h_PR / h_0 = tau litb'" 
response = MsgBox(Msgl3, Style, Title) 
response = MsgBox(msgO, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Style 1, Title) 
itblink = 0 
GoTo 10 
End If 

'Fuel model 1: CSH 
If fuel model name = "CSFI" Then 
f_itb = (1# + f_b) * (tau_litb - tau_lb * tau_tFI) / (temp_fitb) 'eqn(24a) 

Else 

'Fuel model 2: MSH 
loop4 = 0 

Tt4p4 = Tt4 * tau_tFI 'eqn(24b) 

'Call enthalpy(Tt4p4, f_b, h_t4p4)'3/8/05 
Call enthalpy(0, Tt4p4, f_b, h_t4p4, 0, 0) 

f_4p5i = (1# + f_b) * (tau_litb - tau_lb * tau_tFI) / (temp_fitb) 'eqn(24c)initial value for f_4p5i 

4 

If loop4 > 20 Then 
loop4 = 21 
itblink = 0 
GoTo 10 
End If 

'Call enthalpy(Tt4p5, f_4p5i, h_t4p5)'3/8/05 
Call enthalpy(0, Tt4p5, f_4p5i, h_t4p5, 0, 0) 
temp_fitb = eta_itb * h_PR - h_t4p5 
If temp_fitb = 0# Then 
response = MsgBox(Msgl3, Style, Title) 
response = MsgBox(msg0, Style, Title) 

If singlept = False Then jump = MsgBox(MsgJump, Stylel, Title) 
itblink = 0 
GoTo 10 
End If 

f_itb = (1# + f_b) * (h_t4p5 - h_t4p4) / (temp_fitb) 'eqn(24d) 

If Abs(f_b + f itb - f_4p5i) > 0.0001 Then 
f_4p5i = f_b + fitb 
loop4 = loop4 + 1 
GoTo 4 
End If 
End If 

5 

If f_itb = 0# Then 
gamaitb = gamat 
ep_itb = cp_t 
Ritb = R_t 
End If 
'Constants 

cnstl = (gama_itb / R_itb * gc) A (0.5) 

cnst2 = (gama_itb - 1#) / 2# 

cnst3 = (gama_itb + 1#) / (2# * (1# - gama_itb)) 

MFP M8R = M8R * cnstl * (1# + cnst2 * M8R * M8R) A cnst3 'eqn(6r) 

t 

'Component: FAN & LOW-PRESSURE COMPRESSOR PRESSURE RATIOS 
pi_f = (1# + eta_f * (tau_f - 1#)) A cnst8 'eqn(25) 
pi_cL = (1# + eta_cL * (tau_cL - 1#)) A cnst8 'eqn(26) 
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'Component: EXHAUST NOZZLES 
'added at 5/12/04 -1 (Exit Pressure condition) 

'Pt9oP0 = pi_r * pi_d * pi_cL * pi_cH * pi_b * pi_tH * pi itb * pitL * pi_n 
'choked = (0.5 * (gama_itb + 1#)) A (gama_itb / (gama_itb - 1#)) 

'If Pt9oP0 >= choked Then 
' M9 = 1# 

' Pt9oP9 = choked 
' P0oP9 = Pt9oP9 / Pt9oP0 
'Else 

' P0oP9 = 1# 

' Pt9oP9 = Pt9oP0 

' temp_M9 = Pt9oP9 A ((gama_itb - 1#) / gama_itb) 

' If temp_M9 <1# Then 
' response = MsgBox(Msgl4, Style, Title) 

' response = MsgBox(msgO, Style, Title) 

' jump = MsgBox(Msg_jump, Style 1, Title) 

' itblink = 0 
' GoTo 10 
' End If 

' M9 = (2# * (temp_M9 - 1#) / (gama_itb - 1#)) A 0.5 
'End If 

'added at 5/12/04 -2 

Ptl9oP19 = P0oP19 * pi_r * pi_d * pi_f * pi_nf 'eqn(27) 

Pt9oP9 = P0oP9 * pi_r * pi_d * pi_cL * pi_cH * pi_b * pi_tH * pi itb * pitL * pi_n 'eqn(28) 
temp_M9 = Pt9oP9 A ((gama_itb - 1#) / gama_itb) 'eqn(29a) 

If temp_M9 <1# Then 
Msg 14 = " Exit M9 is imaginary! " 

response = MsgBox(Msgl4, Style, Title) 

Msgl4 = " Refer to eqn(27) ==> Pt9/P9 < 1" 
response = MsgBox(Msgl4, Style, Title) 

Msg 14 = "Current values => P0/P9=" & P0oP9 & ", pi_r=" & pi_r & ", pi_d=" & pi_d & ", pi_cL=" & picL 
& ", pi_b=" & pi_b & ", pi_tH=" & pi_tH & ", pi_itb=" & pi_itb & ", pi_tL=" & pi_tH & ", pi_n=" & pi_n 
response = MsgBox(Msgl4, Style, Title) 
response = MsgBox(msg0, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Stylel, Title) 
itblink = 0 
GoTo 10 
End If 

M9 = (2# * (temp_M9 - 1#) / (gama_itb - 1#)) A (0.5) 'eqn(29) 

If M9 >1# Then 
M8= 1# 

Else 

M8 = M9 
End If 

MFP_M8 = M8 * cnstl * (1# + cnst2 * M8 * M8) A cnst3 'eqn(30) 

MFP_M9 = M9 * cnstl * (1# + cnst2 * M9 * M9) A cnst3 'eqn(31) 

A9oA8 = MFP M8 / (MFP M9 * pi_n) 'eqn(32) 

'added at 5/12/04-1 

1 Ptl9oP0 = pi_r * pi_d * pi_f * pi_nf 

1 choked = (0.5 * (gama_c + 1#)) A (gama_c / (gama_c - 1#)) 

1 If Ptl9oP0 >= choked Then 
' M19 = 1# 

1 Ptl9oP19 = choked 
' P0oP19 = Ptl9oP19 / Ptl9oP0 
'Else 

' P0oP19 = 1# 
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' Ptl9oP19 = Ptl9oP0 
' tempM 1 9 = Pt 1 9oP 1 9 A cnst7 
' If gama c < 1# Or temp_M19 <1# Then 
' response = MsgBox(Msgl5, Style, Title) 

' response = MsgBox(msgO, Style, Title) 

' jump = MsgBox(Msg_jump, Style 1, Title) 

' itblink = 0 
' GoTo 10 
' End If 

' M19 = ((temp_M19 - 1#) / cnst5) A 0.5 
' End If 

'added at 5/12/04 -2 

temp_M19 = Ptl9oP19 A cnst7 'eqn(33a) 

If temp_M19 <1# Then 
Msg 14 = " Exit Ml 9 is imaginary! " 

response = MsgBox(Msgl4, Style, Title) 

Msgl4 = " Refer to eqn(31) ==> Pt 1 9/P 1 9 < 1" 

response = MsgBox(Msgl4, Style, Title) 

Msgl4 = "Current values => P0/P19=" & P0oP9 & ", pi_r=" & pi_r & ", pi_d=" & pi_d & ", pi_f=" & pi_f 
& ", pi nf " & pi nl' 

response = MsgBox(Msgl4, Style, Title) 
response = MsgBox(msgO, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Style 1 , Title) 
itblink = 0 
GoTo 10 
End If 

t 

M19 = ((temp_M19 - 1#) / cnst5) A (0.5) 'eqn(33) 

If M19 > 1# Then 
M18 = 1# 

Else 

M18 = M19 
End If 

MFP M18 = M18 * cnst4 * (1# + cnst5 * M18 * M18) A cnst6 'eqn(34) 

MFP M19 = M19 * cnst4 * (1# + cnst5 * M19 * M19) A cnst6 'eqn(35) 

A19oA18 = MFP M18 / (MFP_M19 * pi nf) 'eqn(36) 

t 

'ENGINE BYPASS RATIO 

ebrl = (pi_cLR * pi_cHR * pi_f) / (pi_cL * pi_cH * pi fR) 'eqn(37a) 
ebr2 = (Tt4 * tau rR * tau_fR / (Tt4R * tau r * tau_f)) A (0.5) 'eqn(37b) 
alp = alpR * ebrl * ebr2 * MFP M18 / MFP M18R 'eqn(37) 

t 

'Component: FAN & LOW-PRESSURE COMPRESSOR TEMPERATURE RATIOS 
temp_flpc = tau_cLR - 1# + alp * (tau_fR - 1#) 'eqn(38a) 

If temp_tlpc = 0# Then 

Msgl6 = " Division by zero in eqn(3 8) !!" 
response = MsgBox(Msgl6, Style, Title) 

Msg 16 = "Eqn(38) => Check if 'tau cLR - 1 + alpha*(tau_fR - 1) = O'" 
response = MsgBox(Msgl6, Style, Title) 

Msg 16 = "Current values => tau_cLR= " & taucLR & ", alpha=" & alp & ", tau_fR=" & tau fR 
response = MsgBox(Msgl6, Style, Title) 
response = MsgBox(msg0, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Stylel, Title) 
itblink = 0 
GoTo 10 
End If 
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flpc = (1# - tau_tL) * (1# + f_b + f_itb) / (temp_flpc) 'eqn(38b) 
tau_fN = 1# + (tau_fR - 1#) * eta_mL * tau_litb / tau_r * flpc 'eqn(38)orig 
tau_cL = 1# + (tau_fN - 1#) * (tau_cLR - 1#) / (tau_fR - 1#) 'eqn(39)orig 

I 

'Component: LOW-PRESSURE TURBINE 
tau_tL = 1# - eta_tL * (1# - pitL A ((gama_itb - 1#) / gama_itb)) 'eqn(40)orig 
pi tL = pi tLR * (tau tL / tau tLR) A (0.5) * A4p5oR / A8oR * MFP M8R / MFP M8 'eqn(41)orig 

f 

'Checking: Writing intermediate values of Tt4, Tt4.5, pi_f and picH 
GoTo 39 

Ifpi f - pi fmax >0.001 Then 
overpif = "OVER" 
overpif= overpif+ 1 
Else 

overpi_f = "no" 

End If 

IfpicH - pi_cHmax > 0.001 Then 
overpicH = "OVER" 
overpicFI = overpicFl + 1 
Else 

overpicFl = "no" 

End If 

'Check intermediate values for debugging 
Sheets("test").Cells(liew +1,1) = M0 
Sheets("test").Cells(liew +1,4) = loop2 
Sheets("test").Cells(liew +1,2) = loopif 
Sheets("test").Cells(liew +1,3) = loopich 
Sheets("test").Cells(liew +1,5) = Tt4 
Sheets("test").Cells(liew +1,6) = Tt4p5 
Sheets("test").Cells(liew +1,7)= tau_f 
Sheets("test").Cells(liew +1,8)= tau_cFl 
Sheets("test").Cells(liew +1,9)= pi_f 
Sheets("test").Cells(liew + 1, 10) = picH 
Sheets("test").Cells(liew + 1, 11) =overpi_f 
Sheets("test").Cells(liew + 1, 12) =overpi_cH 
Sheets("test").Cells(liew + 1, 13) =ITBswitch 
Sheets("test").Cells(liew + 1, 14) = reduce 
liew = hew + 1 
39 

If reduce = 2 Or reduce = 4 Then GoTo 20 
tempf = Abs(tau_fN - tau_f) 

If temp f > 0.0001 Then 
tauf = tau_fN 
loop2 = loop2 + 1 
GoTo 2 
End If 

I 

'Skip control limits for throttle hook calculation (PTP) 

'If ctype = 2 Then GoTo 30 '3/12/04 
'If reduce = 2 Or reduce = 4 Then GoTo 20 
If pi fmax > 0# Then 

If pi_f > pi fmax And limit <> 3 And reduce = 1 Then 
pifo = pi_f 
Tt4p5o = Tt4p5 
Tt4p5 = Tt4p5 - 2# 
loopif = loopif + 1 


NASA/TM— 2005-213659 


65 



limit = 3 
GoTo 11 
End If 

err_pif = Abs(pi_f - pi_fmax) 

If limit = 3 And reduce = 1 And err_pif > 0.001 Then 
Dp if = pifo - pi_f 

If Dpi f = 0 Then Dpif = 0.001 'avoid division by zero'3/10/05 
DTt4p5 = Tt4p5o - Tt4p5 
pifo = pi_f 
Tt4p5o = Tt4p5 

DifTt4p5 = (pi_f - pi_fmax) * DTt4p5 / Dpi_f 

If Abs(DifTt4p5) > 100# Then DifTt4p5 = DifTt4p5 / Abs(DifTt4p5) * 100# 
Tt4p5 = Tt4p5 - DifTt4p5 
limit = 3 
GoTo 11 

Elself reduce = 1 Then 
reduce = 2 
GoTo 20 
End If 

If pi_f > pi fmax And limit <> 3 And reduce = 3 Then 
pifo = pi_f 
Tt4o = Tt4 
Tt4 = Tt4 - 2# 
loopif = loopif + 1 
limit = 3 
GoTo 12 
End If 

err_pif = Abs(pi_f - pi_fmax) 

If limit = 3 And reduce = 3 And err_pif > 0.001 Then 
Dp if = pifo - pi_f 

If Dpi f = 0 Then Dpi f = 0.00001 'avoid division by zero 
DTt4 = Tt4o - Tt4 
pifo = pi_f 
Tt4o = Tt4 

DifTt4 = (pi_f - pifmax) * DTt4 / Dp if 

If Abs(DifTt4) > 100# Then DifTt4 = DifTt4 / Abs(DifTt4) * 100# 

Tt4 = Tt4 - DifTt4 
limit = 3 
GoTo 12 

Elself reduce = 3 Then 
reduce = 4 
End If 
End If 
20 

If reduce = 2 Or reduce = 4 Then 
If pi cHmax > 0# Then 
If picH > pi cHmax And limit o 4 Then 
'err_picH = pi_cH - pi cHmax 
If loopich = 0 Then 
pi_cHo = pi_cH 
Tt4o = Tt4 
Tt4 = Tt4 - 2# 

If loopich = 0 Then loopich = 1 
limit = 4 
GoTo 12 
End If 
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End If 
End If 
End If 
30 

Pt3 = P0 * pi_r * pi_d * pi_cL * pi_cH 
Tt3 = TO * taur * taucL * taucH 
'Sheets("SinglePt_FTP").Cells(23, 9) = PRPMhp 
If reduce = 2 Or reduce = 4 Then 
If Pt3max > 0# Then 
If Pt3 > Pt3max And limit <> 5 Then 
Pt3o = Pt3 
Tt4o = Tt4 
Tt4 = Tt4 - 2# 

If loopt3 = 0 Then loopt3 = 1 
limit = 5 
GoTo 12 
End If 

err_Pt3 = Abs(Pt3 - Pt3max) 

If limit = 5 And err_Pt3 > 0.04 Then 
DPt3 = Pt3o - Pt3 
DTt4 = Tt4o - Tt4 
Pt3o = Pt3 
Tt4o = Tt4 

DifTt4 = (Pt3 - Pt3max) * DTt4 / DPt3 

If Abs(DifTt4) > 100# Then DifTt4 = DifTt4 / Abs(DifTt4) * 100# 
Tt4 = Tt4 - DifTt4 
limit = 5 
GoTo 12 
End If 
End If 
End If 

If reduce = 2 Or reduce = 4 Then 
If Tt3max > 0# Then 
If Tt3 > Tt3max And limit <> 6 Then 
Tt3o = Tt3 
Tt4o = Tt4 
Tt4 = Tt4 - 2# 

If loopTt3 = 0 Then loopTt3 = 1 
limit = 6 
GoTo 12 
End If 

err_Tt3 = Abs(Tt3 - Tt3max) 

If limit = 6 And err_Tt3 > 0.1 Then 
DTt3 = Tt3o - Tt3 

If DTt3 = 0 Then DTt3 = 0.00001 'avoid division by zero 
DTt4 = Tt4o - Tt4 
Tt3o = Tt3 
Tt4o = Tt4 

DifTt4 = (Tt3 - Tt3max) * DTt4 / DTt3 

If Abs(DifTt4) > 100# Then DifTt4 = DifTt4 / Abs(DifTt4) * 100# 
Tt4 = Tt4 - DifTt4 
limit = 6 
GoTo 12 
End If 
End If 
End If 


NASA/TM— 2005-213659 


67 



Temp = 100# * (TO * tau r / (TOR * taurR)) A (0.5) 

'PRPMlp = Temp * ((tau_f - 1#) / (tau_fR - 1#)) A (0.5) 'eqn(57) 

PRPMlp = Temp * ((tau_fN - 1#) / (tau_fR - 1#)) A (0.5) ’eqn(57), 3/1 1/05 

PRPMhp = Temp * (tau cL * (tau cH - 1#) / (tau cLR * (tau cHR - 1#))) A (0.5) 'eqn(58) 

If itbinit = 1 Then 'added (see khliew79's note 3/1 1/05) 

If PRPMlpmax > 0# Then 
If PRPMlp > PRPMlpmax And limit <> 7 Then 
PRPMlpo = PRPMlp 
Tt4p5o = Tt4p5 
Tt4p5 = Tt4p5 - 2# 

If loopPRPMlp = 0 Then loopPRPMlp = 1 
limit = 7 

tauf = taufN '3/1 1/05 
GoToll 
End If 

errPRPMlp = Abs(PRPMlp - PRPMlpmax) 

If limit = 7 And err PRPMlp > 0.1 Then 
DPRPMlp = PRPMlpo - PRPMlp 

If DPRPMlp = 0 Then DPRPMlp = 0.001 'avoid division by zero'3/10/05 
DTt4p5 = Tt4p5o - Tt4p5 
PRPMlpo = PRPMlp 
Tt4p5o = Tt4p5 

DifTt4p5 = (PRPMlp - PRPMlpmax) * DTt4p5 / DPRPMlp 
If Abs(DifTt4p5) > 100# Then DifTt4p5 = DifTt4p5 / Abs(DifTt4p5) * 100# 

Tt4p5 = Tt4p5 - DifTt4p5 
limit = 7 
tau_f = tau_fN 
GoTo 11 
End If 
End If 
End If 

If reduce = 2 Or reduce = 3 Or reduce = 4 Then 'modified 3/11/05 
'If reduce = 2 Or reduce = 4 Then 
If PRPMlpmax > 0# Then 
If PRPMlp > PRPMlpmax And limit <> 7 Then 
PRPMlpo = PRPMlp 
Tt4o = Tt4 
Tt4 = Tt4 - 2# 

If loopPRPMlp = 0 Then loopPRPMlp = 1 
limit = 7 
tauf = tau_IN 
GoTo 12 
End If 

err PRPMlp = Abs(PRPMlp - PRPMlpmax) 

If limit = 7 And err PRPMlp > 0.1 Then 
DPRPMlp = PRPMlpo - PRPMlp 

If DPRPMlp = 0 Then DPRPMlp = 0.00001 'avoid division by zero 
DTt4 = Tt4o - Tt4 
PRPMlpo = PRPMlp 
Tt4o = Tt4 

DifTt4 = (PRPMlp - PRPMlpmax) * DTt4 / DPRPMlp 
If Abs(DifTt4) > 100# Then DifTt4 = DifTt4 / Abs(DifTt4) * 100# 

Tt4 = Tt4 - DifTt4 
limit = 7 
tau_f = tau_fN 
GoTo 12 
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End If 
End If 
End If 

If reduce = 2 Or reduce = 4 Then 
If PRPMhpmax > 0# Then 
If PRPMhp > PRPMhpmax And limit <> 8 Then 
PRPMhpo = PRPMhp 
Tt4o = Tt4 
Tt4 = Tt4 - 2# 

If loopPRPMhp = 0 Then loopPRPMhp = 1 
limit = 8 
GoTo 12 
End If 

errPRPMhp = Abs(PRPMhp - PRPMhpmax) 

'If limit = 7 And err PRPMlp >0.1 Then 
If limit = 8 And err PRPMhp >0.1 Then '3/1 1/05 
DPRPMhp = PRPMhpo - PRPMhp 

If DPRPMhp = 0 Then DPRPMhp = 0.00001 'avoid division by zero 
DTt4 = Tt4o - Tt4 
PRPMhpo = PRPMhp 
Tt4o = Tt4 

DifTt4 = (PRPMhp - PRPMhpmax) * DTt4 / DPRPMhp 
If Abs(DifTt4) > 100# Then DifTt4 = DifTt4 / Abs(DifTt4) * 100# 

Tt4 = Tt4 - DifTt4 
limit = 8 
GoTo 12 
End If 
End If 
End If 
itblink = 3 
10 

End Sub 
Sub itbperOQ 
Dim a9, al9, areaO, emf 
Dim fmOl, fm02, fm03 
Title = "ITB Performance Cycle Analysis" 

Style = vbYes 
Stylel =vbYesNo 

t 

'Engine mass flow 

emf = P0 * pi_r * pi_d * pi_cL * pi_cH / (P0R * pi rR * pi_dR * pi cLR * pi cHR) 'eqn(42a) 
mdotO = mdotOR * (1# + alp) / (1# + alpR) * emf * (Tt4R / Tt4) A (0.5) 'eqn(42) 

t 

'Overall engine performance 

fo = (f_b + f_itb) / (1# + alp) 'eqn(43) 

T9oT0 = Tt4p5 / TO * tau_tL / (Pt9oP9) A ((gama_itb - 1#) / gama_itb) 'eqn(44) 

T19oT0 = tau r * tau_f / (Ptl9oP19) A cnst7 'eqn(45) 

a9 = (gama_itb * R_itb * gc * T9oT0 * TO) A (0.5) 'eqn(46) 

al9 = (gama_c * R_c * gc * T19oT0 * TO) A (0.5) 'eqn(47) 

v9 = M9 * a9 'eqn(48) 

vl9 = M19 * al9 'eqn(49) 

v9oa0 = v9 / aO 'eqn(48a) 

vl9oa0 = vl9 / aO 'eqn(49a) 

t 

fmOl = (1# + f_o * (1# + alp)) * v9oa0 + alp * vl9oa0 - (1# + alp) * M0 'eqn(50a) 

fm02 = (1# + fo * (1# + alp)) * R_itb / R_c * T9oT0 / v9oa0 * (1# - P0oP9) / gama_c 'eqn(50b) 
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fm03 = alp * T19oT0 / vl9oa0 * (1# - P0oP19) / gama_c 'eqn(50c) 

FomdO = aO / gc * (fmOl + fm02 + fm03) / (1# + alp) 'eqn(50) 

F = mdotO * FomdO 'eqn(51) 

S = fo / FomdO * units 'eqn(52) 

If Sheets("Input").ComboBoxl. Value = "SI" Then unitchange = 1000000# 

If F < 0# Or S < 0# Then 
wdata = "No" 

Msg 14 = " Thrust or S is negative! " 'modified 3/29/05 
response = MsgBox(Msgl4, Style, Title) 
response = MsgBox(msgO, Style, Title) 

If singlept = False Then jump = MsgBox(Msg_jump, Stylel, Title) 

GoTo 9 
End If 

' added on 7/21/04 -1 
vO = MO * aO 

'fitb = fitb / (1# + f_b) 'then, 'fitb' will be equivalent to 'fitb' in on-design code. This is needed for eta 
computation, (removed 3/8/05) 
fo = f_b + fitb 

corenume = (1# + fo) * v9 * v9 - vO * vO 

fanume = alp * (vl9 * vl9 - vO * vO) 

nume = 0.5 / gc * (corenume + fanume) 

denu = f_b * h_PR * eta b + f itb * h_PR * eta itb 

eta_thermal = nume / denu * 100# 'eqn(54) 

eta_P = (1# + alp) * FomdO * vO / nume * 100# 'eqn(53) 
eta_0 = (1# + alp) * FomdO * vO / denu * 100# 'eqn(54) 

'eta_0 = eta_thermal * eta_P / 100# 

' added on 7/21/04 -2 

areaO = mdotO * TtO A (0.5) / (PtO * MFP MO) 'eqn(56) 

9 

bigit = vbYes 

10 

End Sub 
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